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ABSTRACT 


Transient  cavitation  produced  by  pulsed  ultrasound  in  the  MHz 
range  has  been  investigated.  First,  a  theoretical  background  which 
develops  bubble  motion  in  general  and  focuses  specifically  on  transient 
cavitation  is  presented.  Second,  the  development  of  an  apparatus 
which  is  designed  specifically  for  studying  this  type  of  cavitation  is 
discussed.  Third,  the  use  of  the  apparatus  to  make  measurements  of 
cavitation  thresholds  at  various  frequencies,  pulse  durations,  and  pulse 
repetition  rates  is  discussed.  Finally,  results  of  the  cavitation  thresh¬ 
old  measurements  are  interpreted  using  theories  developed  earlier. 
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I.  INTRODUCTION 


The  term  “cavitation"  was  first  used  in  the  late  nineteenth  cen¬ 
tury  by  Sir  Charles  Parsons  to  describe  the  mechanism  that  caused 
reduced  thrust  in  marine  propellers.  Lord  Rayleigh,  in  1917,  identi¬ 
fied  cavitation  as  the  cause  of  erosion  in  pump  components  that  were  sub¬ 
ject  to  reduced  liquid  pressures.  Cavitation  has  come  to  be  defined  as 
the  rupture  of  a  liquid  or  a  liquid/solid  interface  due  to  the  reduction 
of  liquid  pressure  (Ref.  1J.  Since  the  theoretical  tensile  strength  of 
pure  liquids  is  exceptionally  high,  extremely  large  negative  pressures 
(»  -1000  bars)  should  be  required  to  cause  the  liquid  bonds  to  break. 

Since  rupture  in  real  liquids  can  occur  at  relatively  low  negative 
pressures  (-  -1  bar),  it  has  been  postulated  that  inhomogeneities  exist 
within  the  fluid  which  cause  this  reduced  tensile  strength.  Further¬ 
more.  experimental  results  poiut  to  the  necessity  of  these  micro¬ 
scopic  tnhomogeneities  containing  gas  {Ref.  2J.  Thus,  cavitation  is 
actually  the  formation  of  macroscopic  bubbles  from  microscopic  bub¬ 
bles  (cavitation  nuclei).  Acoustic  cavitation  can  therefore  be  defined  as 
the  activity  of  bubbles  within  a  liquid  that  is  subject  to  an  acoustic 
pressure  field. 

Microcavitation,  or  the  formation  of  micron-sized  bubbles  from 
cavitation  nuclei  due  to  the  influence  of  ultrasonic  acoustic  energy 
(ultrasound),  has  benefits  when  used  In  such  applications  as  ultrasonic 
cleaning  and  etching,  acceleration  of  certain  chemical  reactions,  and. 
medically,  to  destroy  unwanted  cell  formations.  Microcavitation  may. 
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however,  produce  adverse  biological  effects  as  well.  In  particular,  the 
improper  use  of  diagnostic  or  therapeutic  ultrasound  could  possibly 
result  in  the  unintentional  destruction  of  or  damage  to  healthy  cells. 

Investigation  of  microcavitation  phenomena  and  their  possibly 
adverse  biological  effects  has  produced  a  significant  amount  of  infor¬ 
mation  on  cavitation  thresholds.  Robert  Apfel’s  studies  of  pulsed 
ultrasound  cavitation  effects  [Ref.  31  have  summarized  the  possible 
biological  effects  in  single  cells  due  to  excessive  sound  pressure 
amplitudes.  His  work  has  accounted  for  the  effects  of  frequency,  liq¬ 
uid  Inertia,  and  fluid  viscosity  in  the  cavitation  process.  H.  G.  Flynn 
has  done  extensive  numerical  modelling  of  cavitation  phenomena,  also 
taking  into  account  viscosity  and  surface  tension  as  well  as  compress¬ 
ibility  and  heat  conduction  (Ref.  41.  R  A.  Roy.  et  al  (Ref.  51,  and  A.  A. 
Atchley.  et  al.  [Ref.  6).  have  produced  precision  techniques  for  mea¬ 
suring  transient  acoustic  cavitation  thresholds  by  detecting  the  light 
and  sound  emissions  produced  during  violent  bubble  collapse.  It  is 
their  research  that  has  become  the  basis  for  the  research  reported  in 
this  thesis.  Although  some  aspects  of  cavitation  thresholds  have  been 
previously  investigated,  there  has  been  little  work  done  to  determine 
the  time  evolution  of  the  cavitation  event.  The  objectives  of  this 
research  are  to  (1)  develop  an  apparatus  in  which  cavitation  can  be 
induced  and  detected  simultaneously  by  two  methods- acoustic 
scattering  and  sonolumlnescence;  (2)  determine  transient  cavitation 
thresholds;  and  (3)  examine  the  time  evolution  of  the  transient 
cavitation  events. 
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n.  THEORY 


A.  BACKGROUND 

The  behavior  of  bubbles  under  the  influence  of  an  acoustic  pres¬ 
sure  field  is  indeed  complex.  In  this  study,  we  will  examine  the  pres¬ 
sure  amplitudes  of  pulsed  ultrasound  which  will  induce  the  rapid 
growth  and  collapse  of  cavitation  nuclei.  We  will  vary  the  frequency  of 
the  acoustic  signal,  the  number  of  acoustic  cycles  per  pulse,  and  the 
rest  time  between  pulses  in  order  to  more  fully  understand  the  effects 
of  each  parameter  on  the  pressure  amplitude  required  to  cause  the 
nuclei  to  behave  in  this  manner.  In  addition,  we  will  develop  an  appa¬ 
ratus  that  provides  for  both  acoustic  and  optical  observation  of  the 
cavitation  events.  We  must  therefore  be  aware  not  only  of  the  behavior 
of  the  bubbles  when  influenced  by  the  sound  field  but  also  of  their 
propensity  for  scattering  sound  and  producing  visible  light. 

In  this  chapter,  we  will  review  the  above-mentioned  aspects  of 
bubble  behavior.  We  will  first  look  at  the  three  regimes  of  acoustic 
cavitation:  stable  cavitation,  transient  cavitation,  and  gas  body  activa¬ 
tion.  We  will  describe  the  theory  that  describes  bubble  activity  under 
static  equilibrium  conditions,  under  the  influence  of  a  continuous 
sound  pressure  field  in  which  controlled  bubble  oscillations  occur, 
under  conditions  which  result  in  controlled  growth  of  the  bubbles,  and 
under  conditions  which  lead  to  unstable  growth  and  violent  collapse. 
A  theory  will  be  described  that  allows  one  to  predict  the  peak  acoustic 
pressure  which  produces  the  collapse  of  a  gas  bubble.  Finally,  we  will 
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examine  the  theories  which  explain  both  light  emission  (sonolumines- 
cence)  upon  bubble  collapse  and  acoustic  scattering  from  multiple 
bubbles. 

B.  CAVITATION  REGIMES 

It  is  generally  acknowledged  that  two.  and  possibly  three,  cavita¬ 
tion  regimes  exist.  The  focus  of  our  research  will  be  transient  cavita¬ 
tion  produced  by  pulsed  ultrasound.  As  the  name  indicates,  this 
regime  is  transient  in  nature;  bubbles  are  likely  to  exist  for  less  than  a 
few  acoustic  cycles.  The  results  of  the  measurements  to  be  described 
later  indicate  that  effects  of  the  second  regime,  stable  cavitation,  may 
also  be  important.  Stable  cavitation  occurs  when  a  bubble  subjected  to 
a  sound  field  oscillates  about  its  equilibrium  radius.  This  regime  is 
characterized  by  bubbles  that  are  relatively  permanent;  noticeable 
changes  in  the  equilibrium  radius  occur  only  after  many  cycles.  The 
third  regime  (considered  by  many  not  to  be  a  true  cavitation  regime) 
is  known  as  gas  body  activation,  which  involves  the  oscillations  of  gas 
pockets  permanently  embedded  in  cracks,  crevices,  or  within  biologi¬ 
cal  systems.  As  gas  body  activation  does  not  include  actual  formation 
and  growth  of  gas  bubbles,  only  stable  and  transient  cavitation  will  be 
disc  used  further  in  this  chapter. 

C.  BUBBLES  UNDER  QUASI-STATIC  EQUILIBRIUM  CONDITIONS 

The  first  stage  of  any  cavitation  problem  is  nucleation.  or  the  for¬ 
mation  of  a  macroscopically  observable  bubble  from  a  cavitation 
nucleus.  This  formation  usually  Involves  the  loss  of  mechanical 


1 1 


stability  of  the  nucleus.  Following  the  developments  of  Walton  and 
Reynolds  [Ref.  7],  Noltingk  and  Neppiras  [Ref.  8],  and  Blake  [Ref.  9],  a 
brief  examination  will  be  made  of  the  conditions  required  for  a  spher¬ 
ical  bubble  to  lose  its  static  equilibrium  in  a  liquid. 

Given  a  bubble  of  radius  Ro  in  static  equilibrium  with  a  liquid  of 
infinite  extent  at  a  hydrostatic  pressure  Po,  the  pressure  inside  the 
bubble  is  given  by 


p  =  Po+(l? 

/< 2(j\ 

The  term  1^- 


(2.1) 


is  known  as  the  Laplace  or  Young  pressure  and  is  due  to 
the  surface  tension  of  the  liquid,  o  (the  usually  minor  contribution  of 
the  vapor  pressure  has  been  neglected).  Assuming  that  the  liquid 
undergoes  an  incremental  change  in  hydrostatic  pressure  to  a  new 
value  Poo,  the  bubble  radius  will  undergo  an  incremental  change  to 
some  new  radius  R.  Assuming  isothermal  conditions  and  an  ideal  gas. 
we  note 


PiVi  »  P2V2  a  constant 


where 


and 
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p2v2  = 


The  assumption  of  isothermal  conditions  is  justified  because  the  rate 
of  evaporation  and  condensation  of  vapor  at  the  bubble  surface  is  suffi¬ 
cient  to  maintain  this  condition.  After  combining  these  two  expres¬ 
sions  and  cancelling  common  terms,  we  find  that 


Poo  = 


2a 
R  ‘ 


(2.2) 


Figure  2.1.  which  is  a  plot  of  Poo  as  a  function  of  R.  offers  insight 
into  the  stability  of  bubbles.  Assuming  there  is  no  gas  diffusion  into  or 
out  of  bubble.  P»  is  the  value  of  the  liquid  pressure  necessary  to 
balance  the  sum  of  the  internal  pressure  (which  tends  to  make  the 
bubble  expand)  and  the  Laplace  pressure  (which  makes  the  bubble 
collapse).  Poo  is  a  minimum  for  a  bubble  of  radius  R*.  Bubbles  of  radius 
less  than  R*  are  in  stable  equilibrium,  whereas  bubbles  of  radius 
greater  than  R*  are  in  unstable  equilibrium.  To  see  this  more  clearly, 
consider  a  bubble  of  radius  R<R*.  If  R  is  Increased,  P»  decreases 
because  the  decrease  in  the  Laplace  pressure  is  less  than  the 
decrease  in  internal  pressure.  Consequently,  the  external  liquid 
pressure  needed  to  maintain  equilibrium  is  now  less  than  the  actual 
external  pressure.  As  a  result,  the  bubble  shrinks  back  to  R.  In  other 
words,  the  collapsing  pressure  dominates  over  the 
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Figure  2.1 

Graph  of  Pressure  vs.  Bubble  Radius  Showing  Stable 
and  Unstable  Equilibrium  States  of  Bubbles 
Within  a  Liquid  (from  (Ref.  7] 

expanding  pressure.  If  the  radius  had  been  made  smaller.  P„  would 
increase  corresponding  to  a  dominance  of  the  expanding  pressure. 
Again,  the  bubble  returns  to  R.  Therefore,  the  bubble  is  in  a  state  of 
stable  equilibrium— perturbations  in  the  radius  cause  the  bubble  to 
return  to  its  initial  radius.  By  invoking  a  similar  argument  for  bubbles 
of  radius  greater  than  R\  it  is  seen  that  the  equilibrium  is  unstable. 
An  increase  in  radius  causes  P«*  to  Increase  due  to  the  dominance  of 
the  expanding  pressure.  The  bubble  would  increase  in  size 
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indefinitely.  A  decrease  in  radius  causes  the  bubble  to  shrink  to  the 
corresponding  stable  value  of  R.  Because  of  its  role  in  determining  the 
stability  of  bubbles,  R*  is  called  the  critical  radius. 

Now,  consider  a  bubble  of  radius  Ro  <  R*  in  equilibrium  at 
pressure  Po-  According  to  Figure  2.1,  it  is  stable.  If  the  pressure  is 
lowered  to  some  value  Po  -  Pa.  the  bubble  will  grow  until  it  reaches  a 
new  stable  equilibrium  radius  R'  (since  R'  <  R*).  If,  however,  the 
pressure  is  reduced  to  Po  -  Pb.  where  Pb  is  the  Blake  threshold 
pressure,  the  bubble  reaches  the  critical  radius  R*.  Any  further 
reduction  in  pressure  provides  no  equilibrium  value  for  the  bubble 
radius,  and  the  bubble  grows  uncontrollably.  To  solve  for  R\  we  must 
minimize  P«  as  a  function  of  radius. 


Rewriting  equation  (2.2)  slightly. 


2oN 


po+5j  BE-  K* 


2<j 


(2.3) 


then  differentiating  with  respe  t  to  R  and  setting  the  result  equal  to 
zero,  we  find 


(2.4) 


Previously,  we  considered  the  Blake  Pressure  to  be  a  function  of 
initial  bubble  radius.  The  converse  view  is  to  consider  the  initial 
radius  to  be  a  function  of  Blake  pressure.  Accordingly,  the  Blake 
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radius  Rb  is  the  initial  radius  of  a  bubble  that,  when  subjected  to  the 
(acoustic)  Blake  pressure  Pb.  grows  to  a  radius  of  R*.  becomes 
unstable,  and  grows  uncontrollably.  To  find  the  Blake  threshold 
pressure,  we  first  rewrite  R*  (equation  2.4)  in  terms  of  Rb: 


R*  = 


Po +|§)4 

2o  ’ 


1/2 


(2.5) 


Letting 

*■(*>  + H )$ 


(2.6) 


allows  R*  to  be  expressed  more  simply  as 


••■(nr. 


(2.7) 


or 


K  2o 
r?T=  3 


(2.8) 


Setting  Rq  =  Rb.  R  ®  R*.  and  P<»  =  Pq  -  Pb  in  equation  (2.3).  we  find 


_  _  _  2o  2a 

Poo  —  Pq  —  Pb  s  2  “  ^7  • 


(2.9) 


or 


Po  -  Pb  =  - 


4  a 
3  R* 


(2.10) 
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Rewriting  this  expression  using  equations  (2.7)  and  (2.6),  we  find 


8 

Po  ~  pB  =  -  q  o 


3c 


2<p0 + if 


1/2 


(2.11) 


which  can  be  written  in  normalized  form  as 


Pb  = 


Pb 

Po 


=  1  +  qXb 


3Xb  W2 


4(1  -XB)_ 


(2.12) 


where  Xb"^. 

Summarizing  this  discussion,  a  bubble  of  radius  RB  in  stable  equi¬ 
librium  will  lose  its  stability  if  subjected  to  an  acoustic  pressure  of 
amplitude  Pb*  Upon  loss  of  stability,  the  bubble  will  expand  to  many 
times  its  initial  size;  i.e.,  it  will  form  a  macroscopically  observable 
bubble. 


D.  BUBBLES  UNDER  THE  INFLUENCE  OF  A  SOUND  FIELD 

We  begin  the  examination  of  bubble  motion  under  the  influence  of 
a  sound  field  with  a  description  of  bubble  motion  by  Rayleigh  (Ref.  10). 
He  initially  considered  the  collapse  of  an  empty  spherical  cavity 
located  in  an  infinite  mass  of  incompressible  fluid.  The  hydrostatic 
pressure  in  the  liquid  is  assumed  constant.  Rayleigh  observed  that  the 
work  done  by  the  collapsing  bubble  starting  from  rest  had  to  be  equal 
to  the  kinetic  energy  of  the  fluid.  Thus,  the  following  holds  true: 


P«.  (Ro  -  R3)  3  JP  v2(r)  4jir2dr 


(2.13) 
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where  r  is  the  distance  from  the  center  of  the  bubble  to  a  point  in  the 
liquid  that  is  greater  than  the  radius  of  the  bubble  R  and  v(r)  is  the 

4jt 

velocity  of  the  fluid.  Poo  is  the  hydrostatic  pressure  of  the  liquid  and  — 

/  3  \ 

Rq  -  R3  is  the  decrease  in  volume  of  the  spherical  cavity.  The  right- 

^  '  1 

hand  term  is  the  kinetic  energy  of  the  liquid,  j  niv2.  Since  the  fluid  is 
incompressible,  the  change  in  volume  of  the  liquid  is  equal  to  the 
change  in  volume  of  the  bubble  and  the  following  relationship  is  true: 


4jcr2  v(r)  =  constant. 

d4 

Therefore,  4rcr2  v(r)  =  4rcR2  v(R)  or  v2  (r)  =  •py*v2(R).  Letting  v(R)  =  R, 
it  can  be  shown  that 


9  R4  R2 
v2  -  , 


(2.14) 


Substituting  this  result  into  equation  (2.13)  and  performing  the  inte¬ 
gration  on  the  right-hand  side  yields 

Poo  (Ro  "  R^)  3  2tc  p  R3  R2  .  (2.15) 


Solving  for  R2: 

•  2  [*<> 

R2  =  3pP-iR?-1/ 


(2.16) 
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dR 


dR 


Since  R  is  equal  to  ■jj*  ,  it  is  apparent  that  dt  =  — .  As  R  decreases. 


R 


however,  R  is  compelled  to  approach  infinity.  At  this  point,  Rayleigh 
strives  to  avoid  this  obvious  flaw  by  admitting  that  the  initial  assump¬ 
tion  of  zero  internal  pressure  is  invalid.  He  then  assumes  that  the 
cavity  is  filled  with  a  perfect  gas  undergoing  an  isothermal  compres¬ 
sion.  If  this  assumption  is  true,  Boyle’s  law  holds  and  the  pressure 
within  the  bubble  will  increase  until  the  bubble  begins  to  rebound.  If 
fluid  viscosity  is  zero  and  all  other  losses  ignored,  the  bubble  will  con¬ 
tinue  to  oscillate  between  this  minimum  radius  and  the  initial  radius 
Ro.  Again  invoking  the  assumption  that  initial  pressure  inside  the 
bubble  is  zero  or  at  least  some  constant.  Rayleigh  has  derived  the  fol¬ 
lowing  equation  governing  pressure  within  the  liquid,  p.  at  any  point: 


,  R 

R4 

Ro  ' 

1  s3r 

"3r* 

(2.17) 


The  work  of  Plesset,  Poritsky,  and  Noltlngk  and  Neppiras  con¬ 
tributed  to  the  final  formulation  of  an  equation  that  describes  the 
motion  of  a  bubble  under  the  influence  of  an  acoustic  field.  Plesset 
(Ref.  1 1)  extended  Rayleigh's  analysis,  as  did  Poritsky  (Ref.  121  and 
Noltlngk  and  Neppiras  (Ref.  13].  In  1976,  Lauterbom  (Ref.  14)  com¬ 
bined  these  results  with  previously  described  results  to  form  the  final 
characteristic  equation  governing  bubble  motion,  the  RPNNP  equation. 
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Lauterborn  so  named  this  equation  in  honor  of  the  pioneering 
work  done  by  Rayleigh.  Plesset,  Noltingk,  Neppiras,  and  Poritsky.  In 
review,  he  set  the  work  done  by  the  liquid  at  infinity  minus  the  work 
done  by  the  liquid  layer  immediately  adjacent  to  the  bubble  equal  to 
the  kinetic  energy  of  the  liquid  surrounding  the  bubble: 

^  (Pl  -  Poo)  4.t2  dr  =  2jcpR3R2.  (2.18) 

where  Pl  is  the  pressure  in  the  liquid  layer  immediately  surrounding 
the  bubble.  The  reader  should  note  that  P*,  in  this  development  is 
varying  with  time: 


Poo  =  P0  +  P(t) 

where  Po  is  a  static  pressure  about  which  P(t)  varies.  Therefore,  from 
Reference  14,  Pl  can  be  expressed  as 

where  k  Is  the  polytroplc  exponent  of  the  gas  inside  the  bubble,  p  is 
the  shear  viscosity  of  the  liquid,  and  Pv  Is  its  vapor  pressure. 
Differentiating  (2.18)  and  substituting  (2.19)  for  Pl  and  P«,  »  Po  +  P(t). 
we  find  that 


RR  +  |r2 


(2.20) 
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The  following  assumptions  are  Implicit  in  this  comprehensive 
equation  [Ref.  7]: 

1.  The  bubble  remains  spherical  at  all  times; 

2.  Uniform  conditions  exist  within  the  bubble; 

3.  The  acoustic  wavelength,  X,  is  much  greater  than  R; 

4.  No  body  forces  are  present; 

5.  Bulk  viscosity  is  ignored; 

6.  The  liquid  is  incompressible; 

7.  Gas  content  of  the  bubble  is  constant-  no  diffusion  occurs; 

8.  Vapor  pressure,  Pv.  is  constant  during  the  bubble’s  motion. 

Although  many  simplifying  assumptions  have  been  made  in  deriv¬ 
ing  equation  (2.20),  it  is  a  better-than-reasonable  description  of  the 
dynamics  of  a  bubble  subjected  to  an  ultrasonic  sound  field  [Ref.  7].  In 
addition,  the  RPNNP,  or  simply  the  Rayleigh-Plesset,  equation  is  the 
foundation  of  the  theories  discussed  in  the  next  two  sections. 

E.  STABLE  BUBBLE  GROWTH  DUE  TO  AN  ACOUSTIC  PRESSURE 

FIELD 

Stable  cavitation  and  rectified  diffusion,  although  neglected  in  the 
RPNNP  equation,  are  two  related  and  interesting  processes  which  may 
be  Important  in  the  time  evolution  of  transient  cavitation  events.  The 
stable  growth  of  cavitation  nuclei  to  sizes  which  allow  transient  growth 
is  a  likely  mechanism  for  producing  transient  cavitation.  Rather  than 
undertake  a  rigorous  mathematical  development  of  rectified  diffusion 
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as  in  Reference  15,  a  qualitative  description  of  this  phenomenon  will 
be  presented  here. 

The  importance  of  rectified  diffusion,  or  incremental  enlarge¬ 
ment,  was  first  noted  by  Harvey,  et  al.,  in  1944  [Ref.  16].  Blake  [Ref.  9] 
was  the  first  to  attempt  the  derivation  of  a  mathematical  model  in 
1949.  Two  effects,  the  area  effect  and  the  shell  effect,  appear  to  work 
in  concert  to  produce  bubble  growth. 

1.  TlK-AlStt.EffSSt 

When  a  bubble  contracts,  the  concentration  of  gas  within  the 
bubble  increases.  Since  gas  diffuses  from  regions  of  high  concentra¬ 
tion  to  those  of  low  concentration,  the  gas  diffuses  out  of  the  bubble. 
When  the  bubble  expands,  the  concentration  of  gas  decreases,  and  the 
concentration  gradient  reverses,  the  result  being  a  diffusion  of  gas  into 
the  bubble.  The  diffusion  rate  is  proportional  to  the  area  through 
which  gas  may  flow.  Therefore,  more  gas  enters  the  bubble  when  it  is 
expanded  than  leaves  when  it  is  contracted.  The  net  result  is  an  influx 
of  gas  during  one  complete  oscillation  of  the  bubble. 

2.  The 

The  rate  of  diffusion,  as  noted  above,  is  proportional  to  the 
concentration  gradient.  Crum  [Ref.  151  considers  a  spherical  shell  of 
liquid  surrounding  a  bubble.  He  argues  that,  when  the  bubble  contracts, 
the  shell  expands,  with  a  consequent  reduction  in  gas  concentration  near 
the  bubble  surface.  This  increases  the  concentration  gradient  in  the 
outward  direction  over  that  discussed  above,  so  the  diffusion  of  gas  out 
of  the  bubble  is  even  greater.  Conversely,  when  the  bubble  expands, 
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the  shell  volume  decreases,  the  gas  concentration  increases,  and  the 
concentration  gradient  into  the  bubble  increases.  The  overall  result  is 
the  enhancement  of  the  area  effect. 

3.  Threshold  Equation 

There  exist  both  high-  and  low-frequency  thresholds  for 
rectified  diffusion.  Crum  has  formulated  expressions  for  both.  Of 
interest  here  is  the  high-frequency  approximation  [Ref.  15]: 


(  r-,2  2'\ 

PRoCOq 
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2o  Ci 

RoPoo  Co. 


(3+4K) 


(£i)  f 

v  Ro  P»J 

(2.21) 


where 


Pa  *  pressure  threshold  for  rectified  diffusion 


P2  =  p  (D2  Ro/3  Pa* 


b 

2 

(D0 


=  4cop/3  P»  (Damping  term) 
w  4  o 


f  1 
pRo 


3P-  + 


Ro 


(small  amplitude  resonance  frequency 
squared) 


Co  *  equilibrium  concentration  of  gas  in  the  liquid 
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Ci  =  concentration  of  dissolved  gas  in  the  liquid  far  from  the 
bubble 

ji  =  viscosity  of  the  liquid 
©  =  angular  frequency  of  the  acoustic  field 

Ro  =  equilibrium  bubble  radius 
p  =  liquid  density 
Pc  =  ambient  pressure 

*  surface  tension  of  the  liquid 

F.  BUBBLES  DRIVEN  TO  UNSTABLE  GROWTH  AND  VIOLENT 
COLLAPSE 

The  final  step  in  development  of  the  theory  which  defines  bubble 
behavior  In  a  sound  field  is  taken  from  research  by  R  A.  Apfel  and  H.  G. 
Flynn.  In  Reference  3,  Apfel  Indicates  that  the  possibly  adverse  bio¬ 
logical  effects  of  acoustics  are  linked  to  the  growth  and  collapse  of 
submicron  cavitation  nuclei.  A  nucleus  “that  grows  to  5  pm  in  one 
acoustic  cycle  will,  upon  collapse,  deposit  in  the  size  scale  of  one  bio¬ 
logical  cell  approximately  300  million  electron  volts  (MeV)"  of  energy 
(Ref.  31.  This  energy,  in  the  form  of  heat,  results  in  internal  gas  tem¬ 
peratures  in  the  range  of  1000-5000  K.  At  such  high  temperatures, 
effects  such  as  free  radical  formation  and  recombination  can  occur. 
These  effects  can  result  in  the  production  of  visible  light  that  can  be 
detected  optically  (this  will  be  discussed  in  detail  later  in  this 
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chapter).  Apfel’s  work  culminates  in  a  formula  that  predicts  the  peak 

acoustic  pressure  which  produces  a  cavitation  event  (collapse)  given 

the  initial  bubble  radius  Ro  and  maximum  internal  gas  temperature 

5000  K  with  acoustic  frequency,  density,  viscosity,  and  surface  tension 

of  the  liquid  as  parameters.  He  develops  the  relationship  between  the 

maximum  collapse  temperature  T  and  normalized  acoustic  pressure  p 
pA 

•  (Pa  is  peak  acoustic  pressure;  Po  is  ambient  pressure)  by  relating 

p  to  Rmax  (maximum  bubble  radius)  and  Rmax  to  T.  Although  a  rigor¬ 
ous  development  of  these  two  relations  will  not  be  conducted  here. 
Figure  2.2  shows  readily  how  bubble  size  corresponds  to  applied 
acoustic  pressure.  The  reader  is  referred  to  Reference  2  for  a  more 
rigorous  development  of  T  vs.  Rmax  and  Rmax  vs.  p.  It  is  necessary, 
however,  to  examine  Figure  2.2  carefully  to  determine  the  character  of 
the  Rmax  vs.  p  relationship.  The  maximum  bubble  radius  can  be 
thought  of  as  being  composed  of  two  parts:  Ri,  which  occurs  when 
AP,  the  pressure  difference  across  the  bubble.  Is  negative  (outward 
force);  and  Rmax.  which  occurs  when  AP  is  positive,  yet  the  bubble 
continues  to  expand  due  to  outward  liquid  momentum. 

The  time  over  which  maximum  bubble  radius  is  computed  is  also 
noted  in  Figure  2.2.  The  time  t.  during  which  AP  is  negative,  is  the 
algebraic  sum  of  ta.  ti,  and  At. 

t  a  t2  -  ti  -  At  (2.22) 
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and 


At 


=  —  r.,  a  OL.  +  it 
-3  ^  V  AP  +  AP* 


(2.25) 


PA 

where  p  =  p~  and  pb  is  the  normalized  Blake  threshold  pressure  given 
by  equation  (2.12).  The  value  of  x  is  therefore 


”PbJ 
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(2.26) 


The  maximum  bubble  radius  can  now  be  estimated  by 

11/3 


Rmax  2 


i+|(p-D 


(2.27) 


PR3 


From  the  perfect  gas  law,  -j”  «  constant,  and  its  adiabatic  form  PR3y 
=  constant,  it  can  be  shown  that 


r  =  T0  (y  - 1)  . 


(2.28) 


From  equations  (2.26)  and  (2.27),  above,  we  find 

2  -1  [PoTT-TT 

3*0  y  5 


r  =  T0  (y  - 1) 
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(2.29) 


It  is  now  convenient  to  find  frequency  (embedded  In  x).  resulting  in: 
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(2.30) 


1 

.For 

J(p  -  1)  .  J[P-PbJ 

3tcRo  V  p 

L  V  p  4  v  p  J 

3.46 

4li/  2 

2.+  i 

L3  J 

1/3 

+  Ro  V  P  Po  X  ~  0,46 

where 

X  =  pb  +  p-  2  +  V(p  -  1)  (p  -Pb)» 

In  Table  2.1  from  Reference  3  and  Figure  2.3  from  Refer¬ 
ence  4,  the  effects  of  initial  bubble  radius  and  frequency  in  transient 
cavitation  threshold  can  be  seen.  It  is  important  to  note  that  fre¬ 
quency  and  internal  temperature  have  little  effect  on  the  threshold 
pressures  of  smaller  bubbles,  but  become  increasingly  important  as 
initial  bubble  radius  gets  larger.  This  approximate  formulation  is  not 
meant  to  be  a  substitute  for  more  exact  numerical  calculations  such  as 
are  found  in  Reference  4,  but  it  does  give  a  reasonable  idea  of  the 
manner  in  which  the  parameters  of  frequency,  viscosity,  and  initial 
bubble  radius  affect  the  transient  cavitation  threshold  of  a  liquid. 

G.  SONOLUMINESCENCE 

Since  one  of  the  goals  of  this  thesis  is  to  develop  both  optical  and 
acoustical  cavitation  detection  schemes,  it  is  important  to  present 
some  theories  concerning  sonolumlnescence.  Following  the  review  of 
Walton  and  Reynolds  (Ref.  7],  two  primary  theories  as  to  the  origin  of 
sonolumlnescence  have  been  advanced. 
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TABLE  2.1 


PRESSURE  THRESHOLDS  PREDICTED  BY  EQUATION  (2.30) 


Pressure  (bar)  Pressure  (bar) 
T°  E.  f(MHz)  Rq  =  0.1  |im  Rq  =  0.5  fun 


2500 

1.0 

6.4 

2.0 

5000 

6.4 

2.1 

2500 

2.5 

6.4 

2.8 

5000 

6.4 

3.1* 

2500 

5.0 

6.4 

4.2 

5000 

6.4 

5.1 

Figure  2.3 

Threshold  Pressure  Pt  (In  Bars)  for  Transient  Cavitation 
as  a  Function  of  Initial  Radius  Ro  Over  a  Range  of 
Frequency  fX  at  a  Pulse  Width  of  ljis  (Ref.  4] 
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Noltingk  and  Neppiras  [Ref.  13]  proposed  that  sonoluminescence 
is  the  result  of  black-body  radiation  from  the  gas  contained  within  the 
bubble,  which  theoretically  rises  to  extreme  temperatures  as  the  bub¬ 
ble  undergoes  a  rapid  adiabatic  compression  during  its  collapse. 

Griffing  [Ref.  17]  proposed  that  the  formation  of  highly  reactive 
free  radials  in  the  liquid  during  bubble  collapse  were  responsible  for 
rapid  chemical  changes  in  the  fluid,  resulting  in  the  emission  of  light. 

For  the  most  recent  examinations  of  sonoluminescence,  the 
interested  reader  is  referred  to  References  18  through  22.  The  pre¬ 
ponderance  of  evidence  points  to  the  formation  of  free  radicals  that 
subsequently  play  a  role  in  sonoluminescence.  These  free  radicals,  or 
“reactive  molecules  (or  atoms)  with  unpaired  electrons"  [Ref.  23]  may 
be  formed  as  a  result  of  the  dissociation  of  molecules  subjected  to  the 
high  temperatures  produced  by  a  rapid  adiabatic  compression.  Often, 
if  a  free  radical  (with  its  unpaired  electron)  reacts  with  a  molecule 
having  its  full  complement  of  paired  electrons,  the  result  is  another 
free  radical.  This  chain  of  events  will  ts  self-sustaining  until  all  of  the 
suitable  material  has  been  reacted  with,  the  free  radicals  have  com¬ 
bined  with  each  other  to  form  stable  molecules,  or  all  of  the  free 
radicals  have  collided  with  the  wall  of  the  container,  where  they  are 
removed  from  the  process.  According  to  Griffing,  the  emission  of 
light  may  be  due  to  one  of  these  free  radical  reactions  that  happens  to 
be  of  a  chemiluminescent  nature.  This  hypothesis  has  been  confirmed 
by  a  number  of  subsequent  experiments  (Refs.  18-22]. 
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H.  ACOUSTIC  SCATTERING  PROM  BUBBLES 

The  other  detection  method  that  we  hope  to  use  to  investigate 
transient  cavitation  is  the  scattering  of  an  acoustic  signal  (either  the 
source  signal  or  an  alternate  lower-frequency  active  signal)  from 
bubbles  generated  by  such  events.  This  technique  is  expected  to  take 
advantage  of  the  exceptionally  large  scattering  cross-section  of 
resonant  bubbles.  Since  our  measurements  do  not  require  numerical 
solution  of  scattering  parameters,  however,  a  qualitative  explanation 
rather  than  a  rigorous  mathematical  development  will  be  given  here. 
This  theory  is  taken  from  Clay  and  Medwin  (Ref.  24]. 

A  small,  rigid  sphere  whose  density  and  elasticity  are  greater  than 
that  of  the  medium  backscatters  a  very  small  amount  of  the  acoustic 
energy  when  ensonifled  by  acoustic  frequencies  where  the  wavelength 
is  much  greater  than  the  sphere’s  radius.  The  acoustical  cross- 
section  of  the  sphere  is,  in  fact,  much  smaller  than  the  geometrical 
cross-section  because  of  diffraction  of  the  incident  sound  field  around 
the  small,  rigid  body. 

Bubbles,  however,  differ  from  rigid  spheres  in  two  ways:  (1)  the 
specific  acoustic  impedance  (pc)  of  the  gas  within  the  bubble  is  much 
less  than  that  of  the  water,  and  (2)  bubbles  resonate.  Bubbles  are 
extremely  effective  absorbers  and  scatterers  of  sound  at  frequencies 
near  resonance  frequencies.  Scattering  cross-sections  at  resonance 
are  on  the  order  of  400  times  their  geometrical  cross-sections.  (See 
Figure  2.4.) 
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ka  —  -*> 


Figure  2.4 

Ratio  of  Scattering  Cross-Section  to  Geometrical 
Cross-Section  for  a  Bubble  and  a  Rigid  Sphere  (from  (Ref.  24J) 

Transient  cavitation  events  studied  in  this  investigation  will  result 
in  the  generation  of  many  bubbles  of  varying  sizes.  The  effects  of 
multiple  scattering  must  be  considered.  When  scatterers  are  widely 
spaced,  the  scattering  cross-section  per  unit  volume  is  simply  the 
sum  of  the  individual  scattering  cross  sections  of  the  bubbles  within 
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that  volume.  The  bubble  “cloud"  expected  as  a  result  of  the  transient 
cavitation  event  will  be  much  more  densely  packed,  however,  which 
will  cause  the  individual  scattered  fields  of  the  oscillating  bubbles  to 
interact.  The  result  is  that  the  resonance  curve  of  a  single  bubble 
shown  in  Figure  2.4  is  broadened  and  the  scattering  cross-section  of 
an  array  of  bubbles  is  less  than  the  sum  of  the  individual  cross- 
section®.  A  cloud  of  bubbles- bubbles  packed  randomly  in  a  three- 
dimensional  array—  acts  as  a  pressure  release  surface  at  the  face  of  the 
cluster,  thus  reflecting  acoustic  energy  at  this  interface.  It  is  this 
effect  that  is  expected  to  reveal  the  presence  of  acoustic  cavitation. 
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m.  EXPERIMENTAL  APPARATUS 


A.  BACKGROUND 

The  apparatus  with  which  these  cavitation  measurements  will  be 
made  is  based  upon  a  system  developed  by  A.  A.  Atchley,  L.  A.  Frizzell, 
and  R.  E.  Apfel  [Ref.  6].  Our  apparatus  differs  from  theirs  in  three 
major  areas:  (1)  cavitation  detection  can  be  accomplished  by  two 
methods:  optically  (via  sonoluminescence)  and  acoustically,  either 
individually  or  simultaneously;  (2)  this  apparatus  will  rely  on  random 
bubble  nuclei  attaining  threshold  radii  vice  providing  an  artificially 
introduced  uniform  cavitation  nucleus  to  the  sound  field;  and  (3)  this 
apparatus  is  computer-controlled.  The  system  is  made  up  of  two 
mayor  subsystems:  the  fluid  management  system  and  the  electronics 
system.  The  fluid  management  system  will  distill,  filter,  and  degas 
ordinary  tap  water  for  use  in  the  cavitation  enclosure.  The  electronics 
system  will  provide  the  capability  to  induce,  detect,  and  record  tran¬ 
sient  cavitation  events,  all  under  computer  control. 

B.  CAVITATION  ENCLOSURE  AND  LIGHT-TIGHT  ENCLOSURE 

The  cavitation  enclosure  (Figure  3.1)  is  a  25  cm  by  15  cm  by  15 
cm  box  made  of  1.2  cm  thick  plexiglass.  Plexiglass  was  chosen  as  the 
primary  structural  material  because  it  is  sufficiently  strong  to  support 
Installed  transducers  and  fluid  inlet  and  outlet  connections,  yet  it 
allows  the  optical  detection  apparatus  to  “see"  what  occurs  in  the 
insonified  region  of  the  system.  One  end  of  the  chamber  holds  a 
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Figure  3.1 

Cavitation  Enclosure 


Figure  3.2 

Light-Tight  Enclosure  and  Photomultiplier  Tube  Assembly 


Harisonic  Labs  focused  transducer  that  supplies  the  cavitation  genera¬ 
tion  signal.  A  wedge  and  two  vertical  walls  (also  plexiglass)  have  been 
installed  at  the  opposite  end  of  the  enclosure  to  scatter  the  incident 
sound  and  thus  minimize  reflections  that  would  interfere  at  the  focal 
region.  The  source  transducers  (rated  at  1.0  MHz  and  2.25  MHz)  have 
focal  lengths  of  7.6  cm.  Two  additional  walls  with  holes  cut  to  allow 
passage  of  the  focused  sound  beam  are  installed  on  either  side  of  the 
focal  region,  again  to  minimize  interference  from  reflections.  The 
mounting  for  the  source  transducer  allows  for  rapid  changing  of  ele¬ 
ments  so  that  the  other  source  frequency  can  be  examined.  Two 
Panametrics  V301  transducers  are  permanently  installed  in  the  bot¬ 
tom  and  one  side  of  the  enclosure  in  such  a  manner  that  their  beam 
pattern  intersects  the  focal  region  of  the  source  transducer.  The  top 
of  the  enclosure  is  removable  to  allow  access  to  the  interior.  The 
cavitation  enclosure  is  placed  in  a  light-tight  box  that  has  connections 
for  fluid  inlet  and  outlet  as  well  as  BNC  connections  for  source  and 
receive  transducer  signals  (Figure  3.2).  The  photomultiplier  tube 
(PMT)  assembly  is  mounted  on  one  side  of  the  light-tight  enclosure  so 
that  it  also  can  “see’*  the  focal  region  for  sonoluminescent  indications 
of  cavitation.  The  photomultiplier  tube's  sensitive  circuitry  is  pro¬ 
tected  by  a  mechanical  shutter  that  can  be  lowered  to  block  light  and 
by  an  Interlock  switch  in  the  lid  of  the  enclosure  that  Interrupts  the 
high-voltage  power  supply  to  the  PMT  when  the  lid  is  removed. 
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C.  FLUID  MANAGEMENT  SYSTEM 

The  fluid  management  system  (Figure  3.3)  consists  of  a  four-ele¬ 
ment  deionization/filtration  system,  a  bag-type  0.2  pm  filter,  a  small 
variable-speed  electric  centrifugal  pump,  an  ultraviolet  sterilizer,  and 
an  Erlenmeyer  flask  with  connections  for  system  evacuation  and  argon 
infusion.  The  Erlenmeyer  flask  is  also  used  as  a  fluid  reservoir  and 
expansion  tank.  The  filtration/deionization  system  produces  the  equiv¬ 
alent  of  quintuply  distilled  water  and  removes  all  ionized  minerals, 
silicon,  and  free  carbon  dioxide  down  to  four  parts  per  billion.  The 
bag  filter  ensures  that  no  particles  larger  than  0.2  pm  are  introduced 
into  the  system.  These  two  elements  are  used  only  when  the  system 
is  being  filled.  When  the  system  is  fully  charged,  the  water  is  circu¬ 
lated  through  the  cavitation  enclosure  via  a  closed  loop  consisting  of 
the  pump,  the  U/V  sterilizer,  the  enclosure,  and  the  expansion  reser¬ 
voir.  Argon  is  bubbled  into  the  fluid  to  help  enhance  the  sonolumines- 
cent  reaction.  The  system  is  simultaneously  evacuated  (15  inches  Hg 
vacuum)  in  order  to  remove  dissolved  gases.  The  evacuation  is  contin¬ 
ued  until  visible  air  bubbles  are  removed  (usually  about  one  hour).  The 
circulation  system  is  then  secured  and  the  sample  is  returned  to 
atmospheric  pressure,  where  it  remains  quiescent  for  one-half  hour 
before  data  accumulation  begins.  Alter  about  an  hour  of  data  accumu¬ 
lation  (usually  20  to  40  data  points),  the  system  must  be  drained,  refilled, 
and  reprocessed.  Refilling  the  system  is  necessary  because  the  fluid 
becomes  “tired"  after  multiple  cavitation  events,  causing  higher-than- 
normal  thresholds  (possibly  due  to  fewer  remaining  cavitation  nuclei). 
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Figure  3,3 

Fluid  Management  System 


D.  CAVITATION  GENERATION  AND  DETECTION  SYSTEM 

The  electronics  system  consists  of  the  cavitation  generation  sys¬ 
tem  and  the  cavitation  detection  system  (Figure  3.4).  The  heart  of  the 
cavitation  generation  system  is  the  Hewlett-Packard  9000  Series  300 
personal  computer,  which  controls  the  output  frequency,  pulse  width 
(duty  cycle),  pulse  repetition  rate,  and  acoustic  pressure  amplitude  of 
the  source  transducer  via  the  Hewlett-Packard  3314A  function  gener¬ 
ator.  The  output  of  the  function  generator  is  boosted  by  an  Amplifier 
Research  100-A15  Power  Amplifier  and  is  then  applied  to  the  source 
transducer.  A  printout  of  the  control  program  is  contained  in 
Appendix  A 

The  cavitation  detection  system  consists  of  two  subsystems:  the 
acoustic  detection  system  and  the  photomultiplier  tube  (PMT)  assem¬ 
bly.  The  acoustic  detection  system  is  made  up  of  the  two  Paname tries 
V301  transducers  installed  in  the  cavitation  enclosure  and  the  associ¬ 
ated  elements  for  either  passive  or  active  (via  scattering  of  a  cw  signal) 
detection  of  cavitation  events.  In  the  passive  mode,  the  output  of  one 
of  the  detection  transducers  is  boosted  by  a  small  RF  amplifier  and  is 
displayed  on  the  Tektronix  2245A  oscilloscope.  The  trace  of  this 
signal  on  the  oscilloscope  consists  primarily  of  a  bundle  of  acoustic 
pulses  that  have  either  been  diffracted  or  scattered  to  the  detection 
transducer  face.  The  bundle  is  delayed  from  the  source  pulse  by 
approximately  100  microseconds,  which  corresponds  to  the  time  the 
pulse  would  take  to  make  the  15  cm  journey  from  the  source  to  the 
focal  region  and  then  to  the  detection  transducer  face.  The 
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system  operator  must  monitor  the  signal  in  the  vicinity  of  this  time 
delay  for  a  transient  that  would  indicate  an  increase  in  the  scattered 
source  signal,  thus  signifying  the  occurrence  of  a  cavitation  event.  The 
operator  must  then  manually  interrupt  the  computer  ramp  function 
for  the  computer  to  calculate  and  display  the  threshold  pressure. 

In  the  active  mode,  the  second  Panametrics  transducer  is  driven 
at  a  continuous  frequency  in  the  100  kHz-500  kHz  range.  The  result¬ 
ing  active  signal  is  received  by  the  detection  transducer  along  with  the 
source  signal.  The  output  of  the  detection  transducer  is  passed 
through  a  low  pass  filter  to  remove  the  cavitation  generation  signal 
while  passing  the  active  detection  signal.  A  sharp  increase  in  the 
amplitude  of  the  received  signal  indicates  scattering  from  the  bubble 
cloud  that  is  formed  by  a  cavitation  event.  The  photomultiplier 
assembly  and  associated  circuitry  provide  optical  detection  of  the 
sonoluminescence  generated  by  a  transient  cavitation  event.  The 
pulse  from  the  PMT  is  pre-amplified  and  sent  to  an  ORTEC  9302 
pulse  height  discriminator,  which  determines  the  sensitivity  of  the 
mechanism.  The  discriminated  signal  is  then  sent  through  a  Canberra 
2071 A  counter.  Because  of  background  light  noise  within  the  light¬ 
tight  enclosure,  there  is  always  some  pulsed  output  from  the  PMT. 
Since  transient  events  nearly  always  occur  in  groups,  events  will  likely 
produce  numbers  of  pulses  significantly  higher  than  background.  At 
the  time  of  this  writing,  attempts  to  detect  cavitation  via  this  method 
have  proved  unsuccessful. 


E.  TRANSDUCER  CALIBRATION  APPARATUS 


In  addition  to  the  apparatus  used  for  the  generation  and  detection 
of  cavitation,  the  apparatus  used  for  the  source  transducer  calibration 
is  sufficiently  involved  to  warrant  further  description.  Calibration 
involves  measuring  the  transverse  and  axial  beam  pattern  of  the  source 
transducer  and  determining  the  transducer’s  total  acoustic  power 
output  by  force  target  deflection. 

The  axial  and  transverse  beam  patterns  of  both  the  source  and 
detection  transducers  were  found  using  a  small  FVDF  probe 
hydrophone  designed  and  built  specifically  for  this  use.  The 
construction  and  calibration  of  this  hydrophone  were  accomplished  to 
fulfill  requirements  for  a  final  project  in  a  transducer  theory  and 
design  course.  Pertinent  portions  of  the  report  submitted  on  the 
project  are  contained  on  Appendix  B.  The  PVDF  hydrophone  was 
mounted  in  a  micrometer  Jig  that  allowed  precise  measurement  of  its 
position  relative  to  the  transducer  face  (Figure  3.5). 

The  force  target  deflection  measurement  was  used  to  determine 
peak  output  frequencies  and  total  acoustic  power  output  of  the  source 
transducers  (Figure  3.6).  A  description  of  the  theory  behind  the  force 
target  procedure  follows  in  Chapter  IV.  In  this  procedure,  the  source 
transducer  is  driven  by  the  function  generator  via  the  power  amplifier 
at  the  same  gain  setting  as  will  be  used  in  the  threshold  experiment. 
The  target  is  a  small  plexiglass  plate,  suspended  at  a  45-degree  angle 
by  four  pieces  of  two-pound  test  monofilament  fishing  line  of  equal 
length  (approximately  .04  mm  diameter,  225  pgrams /cm  linear 
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FVDF  Hydrophone  Determination  of  Source  Transducer  Beam 

Pattern 


Figure  3.6 


Equipment  Arrangement  for  Force  Target  Deflection 

Measurements 
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density).  The  plexiglass  target  Is  rectangular  in  shape,  of  dimensions 
3.8  x  5.1  x  0.7  cm.  It  has  been  carefully  weighted  using  lead  strips  to 
achieve  a  weight  of  1  g  in  water  (14.7  g  in  air).  A  cavity  approximately 
0.3  cm  deep  is  machined  into  the  target  and  covered  by  a  clear,  thin, 
stiff  plastic  to  form  a  pressure-release  surface.  This  ensures  that  the 
radiation  force  is  almost  entirely  transformed  into  a  deflecting  force. 
Essentially  no  acoustic  energy  is  reflected  back  to  the  driving 
transducer,  but  rather  toward  the  bottom  of  the  test  tank.  The 
deflection  was  measured  by  a  graduated  telescope  mounted  to  the 
micrometer  jig  described  above. 
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IV.  EXPERIMENTAL  PROCEDURE 


A.  BACKGROUND 

Prior  to  data  taking,  it  is  necessary  to  determine  the  performance 
characteristics  of  the  transducers  involved.  This  involves  reciprocity 
calibration  of  the  two  Panametrics  V301  transducers  used  for  the 
acoustic  detection  of  cavitation  events.  Also,  once  the  source 
transducers’  true  peak  output  frequencies  are  determined,  it  is 
necessary  to  determine  their  axial  and  transverse  beam  patterns. 
These  transducers  are  used  to  generate  the  ultrasonic  field  that 
induces  cavitation.  The  axial  beam  patterns  are  used  to  verify  the 
transducers’  focal  lengths  while  the  transverse  beam  patterns  are 
integrated  to  yield  the  effective  area,  A^fj,  which  relates  peak  intensity 
and  total  transducer  output  power. 

Having  determined  the  beam  patterns,  it  is  necessary  to  perform 
a  radiation  force-target  deflection  measurement  of  the  Harisonic 
transducers  in  order  to  calculate  transducer  output  power.  Having 
experimentally  determined  A^ff  and  total  power,  the  output  peak 
intensity  and  consequently  the  output  pressure  characteristics  can  be 
accurately  determined. 

Upon  completion  of  transducer  source  strength  calibration,  cavi¬ 
tation  samples  are  prepared  and  threshold  measurements  are  made  in 
order  to  establish  a  data  base. 


B.  RECIPROCITY  CALIBRATION 


The  technique  employed  to  find  the  sensitivities  of  the  Panamet- 
rics  V301  transducers  is  identical  to  that  described  in  Reference  25. 
See  Figure  4.1.  The  distance  between  transducers  is  20.2  cm.  This 
distance  is  well  into  the  far  field  of  the  transmitter,  which  is  found 
using  the  following  approximation: 


Tmin  1  L 

“~4r 


(4.1) 


where 

X  =  acoustic  wavelength  1.5  mm  (assuming  c  =  1481  m/s  for 
fresh  water  and  f  a  986  kHz), 

L  a  diameter  of  V301  transducer  *  3.2  cm. 


Train  =  minimum  distance  to  the  acoustic  far  field  (in  cm). 

By  employing  the  appropriate  equations  established  in  Reference  25. 
graphs  of  the  transducers’  sensitivities  can  be  constructed.  See  Figure 
4.2.  The  V301  transducers  with  serial  numbers  596  and  598  were  used. 

C.  AXIAL  BEAM  P  iTTERN  MEASUREMENT  OP  SOURCE 
TRANSDUCERS 

Using  the  polyvinylidene  fluoride  (PVDF)  hydrophone  discussed 
above  (see  Appendix  B)  and  a  precision  measuring  Jig  (see  Figure  3.5), 
the  axial  beam  patterns  of  the  two  Harlsonic  source  transducers  are 
determined.  A  CW  sinusoid  at  the  transducer’s  rated  operating 
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frequency  is  generated  by  the  HP  3314A  function  generator  and 
passed  through  the  Amplifier  Research  model  100A15  power 
amplifier  and  finally  to  the  test  transducer.  The  amplifier  gain  is  set  at 
200  for  a  constant  15  VAC  (peak-to-peak)  output.  The  PVDF 
hydrophone  is  inserted  into  the  mounting  bracket  attached  to  the 
measuring  jig.  The  active  element  of  the  hydrophone  is  oriented  such 
that  the  hydrophone  face  moves  along  the  transducer’s  longitudinal 
axis  as  the  micrometer  knob  is  rotated.  By  recording  the  received 
voltage  from  the  hydrophone  and  the  position  from  the  measuring  jig’s 
micrometer  knob,  an  axial  beam  pattern  can  be  determined  and  the 
transducer  focal  length  verified.  The  Harisonic  transducers  used  in 
this  series  of  experiments  were  found  to  have  focal  lengths  of  7.29  and 
7.37  cm.  There  is  not  a  large  enough  difference  between  these  values 
and  the  rated  focal  length  of  7.62  cm  (3  in.)  to  warrant  redesign  of  the 
cavitation  tank  or  relocation  of  the  transducers  used  for  acoustic 
detection  of  cavitation. 

D.  TRANSVERSE  BEAM  PATTERN  MEASUREMENT  OP  SOURCE 

TRANSDUCERS 

1.  Theorr 

It  is  necessary  to  relate  peak  intensity  and  total  acoustic 
power  in  order  to  establish  a  relationship  between  transducer  input 
voltage  and  transducer  output  pressure.  In  order  to  establish  this 
relationship  .  we  define  an  effective  beam  area  Ae#  by 


where  Ipeak  =  peak  intensity  (W /m2)  and  Ptot  =  total  output  power 
(W).  Since  I  a  p2  a  V2,  I  =  K’V2  where  K'  is  a  constant;  p  is  acoustic 
pressure  and  V  is  the  corresponding  hydrophone  ouput  voltage.  This 
voltage,  a  function  of  distance  from  the  center  of  the  transducer  face, 
is  determined  by  a  transverse  scan  of  the  focal  region.  At  the  focal 
point,  Ipeak  =  K'Vpeak-  Substituting  this  expression  into  equation 
(4.2)  gives 


Aeff  = 


Ptot 


PEAK 


(4.3) 


It  is  known  that 


Ptot  =  J  IdS;  dS  *  incremental  area  normal  to  the  beam. 


The  integral  is  taken  over  the  entire  cross-sectional  area  of  the  beam 
and  can  be  approximated  (in  polar  coordinates)  by 


Ptot 


■;r 


K'V2rdrd8  =  2itK 


V^rdr 


(4.4) 


where  r  =  radial  distance  from  transducer  axis  and  a  is  some  radial 
distance  within  which  is  contained  practically  all  of  the  beam  energy. 
In  practice,  a  is  chosen  equal  to  the  radius  of  the  transducer  face. 
Therefore, 


2xK*J  V^rdr  2rtJ  V^rdr 


PEAK 


w 


(4.5) 


PEAK 
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Finally,  since  Aeff  can  be  determined  from  the  transverse  beam  profile, 
it  is  necessary  only  to  calculate  Ptot  (as  described  in  a  subsequent 
section)  in  order  to  calculate  the  peak  intensity  since  Ipeak  = 
FhOTMeff- 

Further,  since  peak  pressure,  Ppeak.  Is  related  to  peak 
intensity  by  the  following  relationship  (assuming  either  plane  or 
spherical  waves), 

Ppeak  =  (pocIpeak)  1/2.  (4.6) 

Ppeak  can  be  rewritten  in  terms  of  and  Ptot-' 

Ppeak  =  PtotJ1/2-  (4.7) 

2.  Prpcfidurc 

The  procedure  used  for  measuring  transverse  beam  patterns 
is  very  similar  to  that  described  in  the  previous  section.  The  Jig  is 
positioned  in  a  manner  such  that  the  hydrophone  travels  in  a  line 
perpendicular  to  the  longitudinal  axis  at  a  distance  from  the 
transducer  face  equal  to  the  focal  length.  Accurate  position  measure¬ 
ments  (±10  pm)  can  be  made  by  recording  the  micrometer  knob 
readings.  Concurrently,  hydrophone  output  voltages  are  recorded 
from  the  oscilloscope  display,  resulting  in  the  transverse  beam  pat¬ 
terns  of  Figures  4.3  and  4.4.  Plots  of  V^r  vs.  r  are  included  as  Figures 
4.5  and  4,6.  Only  the  right  half  of  each  beam  pattern  was  used  to 
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Output  of  PVDF  Hydrophone  (mV) 


Figure  4.3 

Transverse  Beam  Patters  of  1.0  MHz  Transducer 


Olstance  from  Center  of  Transducer  Face  (cm) 


Figure  4.4 

Transverse  Beam  Pattern  of  2.25  MHz  Source  Transducer 
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Distance  from  Center  of  Transducer  Face  (cm) 
Figure  4.5 

Var  vs.  r  (integrand  of  Equation  (4.4)); 
1.0  MHz  Source  Transducer 


Distance  from  Center  of  Transducer  Face  (cm> 


Figure  4.6 

V?r  va.  r  (Integrand  of  Equation  (4.4)); 
2.25  MHz  Source  Transducer 
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be  symmetrical.  This  is  not  exactly  true;  but  no  appreciable  error  is 
introduced  by  making  this  assumption.  Integrations  of  these  plots  are 
performed  by  using  a  Macintosh  Plus  personal  computer  and  DBGPLOT 
software.  The  results  are  plotted  in  Figures  4.7  and  4.8. 

E.  DETERMINATION  OP  TRANSDUCER  PEAK  OUTPUT  FREQUENCY 

In  order  to  optimize  the  output  from  the  source  transducers,  it  is 
necessary  to  know  their  peak  output  frequency.  Using  the  deflection 
apparatus  illustrated  in  Figure  3.6.  the  rectangular  target  is  suspended 
in  the  focal  region  of  a  source  transducer.  A  sinusoidal  signal  at  the 
rated  operating  frequency  is  applied  to  the  transducer  via  the  power 
amplifier.  A  constant  input  voltage  of  15  volts  (peak  to  peak)  is  main¬ 
tained  (monitored  on  the  oscilloscope)  while  frequency  is  varied  above 
and  below  the  manufacturer's  quoted  peak  output  values.  Target 
deflection  is  measured  for  the  various  test  input  frequencies  by  noting 
the  change  in  target  position  with  a  graduated  telescope  mounted  on 
the  measuring  Jig  used  previously.  Plots  of  target  deflection  versus 
transducer  frequency  are  reproduced  as  Figures  4.9  and  4. 10.  The 
transducers  specified  as  1  MHz  and  2.25  MHz  peak  output  frequencies 
were  found  to  have  true  peak  output  frequencies  of  986  kHz  and  2.17 
MHz.  respectively. 


•  eA  i-OtX)  jp  j  ZA  / 
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0.0  0.5  1.0  15  2.0  2.5 

Distance  from  Center  of  Transducer  Face  (cm) 


Figure  4.7 

Integral  of  V^rdr  for  1.0  MH*  Transducer 


Distance  from  Center  of  Transducer  Face  (cm) 

Figure  4.8 

Integral  of  V2rdr  for  2.25  MHz  Transducer 
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Transducer  Operating  Frequency  (MHz) 

Figure  4.9 

Target  Deflection  versus  Frequency  for  Determination 
of  Peak  Operating  Frequency  of  1.0  MHz  Source  Transducer 


Figure  4.10 

Target  Deflection  versus  Frequency  fo^  Determination 
of  Peak  Operating  Frequency  of  2.25  MHz  Source  Transducer 
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F.  RADIATION  FORCE-TARGET  DEFLECTION  MEASUREMENT 


1.  Theory 

The  method  used  to  determine  the  peak  output  pressure  of 
the  Harisonic  transducers  is  a  variation  of  the  method  described  by 
Dunn,  et  al.  [Ref.  26].  Rather  than  using  a  sphere,  as  in  Dunn’s  work, 
the  rectangular  plexiglass  target  is  employed. 

A  rectangular  target  is  used  rather  than  a  sphere  because,  as  a 
sphere  deflects,  it  tends  to  move  out  of  the  sound  field,  while  a  rela¬ 
tively  large  plate  is  less  susceptible  to  this  effect.  In  addition,  there  is 
no  need  to  calculate  the  form  function. 

Referring  to  Figure  4.11  (reproduced  from  Reference  26  with 
some  clarifying  remarks),  it  is  a  relatively  simple  matter  to  construct 
the  appropriate  free-body  diagram  and  the  equations  governing  the 
static  equilibrium: 

IFX  =  -Fr  +  Tcos0  =  0. 

£Fy  =  -mg  +  Tsin0  =  0. 


The  above  equations  imply  that 


T  Fr  nig 
“  cos0  “  sln0’ 

Therefore.  Fr  =  However.  tan0  =  If  L»h.  tan0 

Therefore. 


(4.8) 

L/d. 


~  tan0  "  “  L  ’ 


(4.9) 


$6 


Tramductr 


Figure  4.11  — 

Free  Body  Diagram  of  Force  Target  Deflection 

From  References  25  and  27,  it  can  be  shown  that  the  excess  pressure 
pe  -  pc2s,  where  p  =  density  of  the  medium,  s  =  condensation,  and  c  a 

phase  speed  of  the  acoustic  waves.  For  a  high-intensity  sound  field, 
variations  in  density  are  significant  so  p  s  poll  ♦  s).  By  substitution. 


pe  =  pc2s  =  po(l  +  s)c2s  =  poc2s  -4-  poc2s2.  If  s  =  socoscot,  the  time- 
averaged  excess  pressure,  <pe>t*  can  be  rewritten  as  follows: 

2 

<Pe>t  =  POC2S0  <COSO)t>  +  PoC2Sq  <cos2cot>. 

The  first  term  on  the  right-hand  side  reduces  to  zero,  but  the  second 
term  (or  radiation  pressure)  can  be  expanded  by  employing  the 
appropriate  trigonometric  identity. 

pOC2S0  <COS2t*:i>  =  POC2Sq  •  2  [1  +  <COS2C0t>J 

2  2 
P0C2s0  Poc2s0 

«  — ^ — +  — 2 —  <cos2©t>. 

2  POC2Sq  j 

This  reduces  finally  to  poc2sq  <cos2cat>  *  — g —  •  However,  pe  »  — 

where  I  =  intensity.  Recalling  that  power.  Ptot  =  I  •  A  where  A  =  area, 

it  can  be  seen  that  the  relationship  Ptot  =  <Pe>tcA  is  also  true. 

Setting  <pc>t  •  A  **  Fr,  Ptot  *  cA  *  Fr  •  c  where  c  »  1481  m/s  since 

these  tests  are  conducted  in  fresh  water.  By  substitution, 

Ptot  a  •  d  (4.10) 

mcc 

where  — f  -  is  a  constant  and  d  is  the  measured  target  deflection. 

By  substituting  the  appropriate  values  for  m,  g,  c.  and  L. 


This  equation  is  used  in  the  radiation  force-target  deflection  mea¬ 
surements.  Upon  completion  of  such  measurements,  all  ingredients 
necessary  to  characterize  transducer  output  pressure  will  have  been 
determined. 

2.  Procedure 

A  CW  sinusoid  (at  the  transducer’s  peak  output  frequency) 
generated  by  the  HP  3314A  function  generator  is  sent  through  the 
Amplifier  Research  model  100A15  power  amplifier  and  finally  to  the 
Harlsonic  transducer  undergoing  test.  The  deflection  of  the  plexiglass 
target  is  measured  as  the  transducer  input  voltage  (as  monitored  on 
the  Tektronix  2445A  oscilloscope)  is  varied.  Plots  of  target  deflection 
vs.  input  voltage  squared  are  included  as  Figures  4,12  and  4.13.  The 
slopes  of  these  plots  are  proportional  to  output  power  for  a  given  input 
voltage.  Note  that  the  “mirror’*  Images  of  the  recorded  data  points  are 
plotted  as  well.  This  ensures  that  a  line  determined  by  linear 
regression  will  pass  through  the  origin,  resulting  in  zero  power  output 
for  zero  input  voltage,  as  is  dictated  by  common  sense. 

Having  determined  all  of  the  terms  necessaiy  to  characterize 
output  pressure,  the  pressure  as  a  function  of  input  voltage  can  be 
plotted.  See  Figures  4.14  and  4.15.  Again,  the  “mirror"  images  are 
plotted  as  well  in  order  to  drive  the  linear  regression  through  the 
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Figure  4.14 

Output  Pressure  Amplitude  vs.  Input  Voltage 
for  1.0  MHz  Source  Transducer 


Transducer  Input  Voltage  (Volts) 

Figure  4.15 

Output  Pressure  Amplitude  vs.  Input  Voltage 
for  2.25  MHz  Source  Transducer 
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origin.  Note  that,  in  Figure  4.15,  the  relatively  large  deviations  from 
the  linear  regression  at  lower  voltages  are  of  no  concern  since  the 
transducer  is  not  operated  in  this  region.  The  relationships  between 
the  peak-to-peak  transducer  input  voltage  and  the  peak  output 
pressure  are  included  below. 

For  the  1  MHz  transducer: 

Ppeak  (bar)  =  0.0522  Vin  (mV). 

For  the  2.25  MHz  transducer: 

Ppeak  (bar)  =  0.0403  Vin  (mV). 


G.  DETERMINATION  OP  TRANSIENT  CAVITATION  THRESHOLDS 

All  threshold  data  were  taken  using  the  passive  acoustic  (scatter¬ 
ing  of  the  source  signal  from  the  bubble  “cloud*  produced  by  a  tran¬ 
sient  event)  detection  scheme.  This  method  was  the  first  to  show 
acceptable  reliability.  As  the  data  required  was  substantial  and  excep¬ 
tionally  time-consuming  to  take,  we  concentrated  on  data  accumula¬ 
tion  via  this  method  rather  than  sacrifice  more  time  refining  the  other 
two  detection  methods.  The  HP  9000  series  300  computer  was  pro¬ 
grammed  to  allow  user  input  of  desired  frequency,  pulse  repetition 
frequency  (PRF),  pulse  duration,  and  pressure  amplitude  sweep  para¬ 
meters  (via  input  of  function  generator  starting  amplitude,  voltage 
amplitude  increase  per  step,  and  total  number  of  amplitude  steps). 
Data  was  obtained  at  986  kHz  and  2.17  MHz  using  200  Hz,  500  Hz, 
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and  1000  Hz  PRF.  Pulse  durations  of  1.5,  2,  5,  10,  15,  20,  25,  30,  35, 
and  40  psec  were  examined.  All  voltage  ramps  were  started  at  300- 
400  mV  peak  input  voltage  to  the  power  amplifier,  which  results  in  an 
initial  peak  output  pressure  from  the  source  transducers  of  approxi¬ 
mately  3-4  bar.  Incrementally  increasing  pressure  steps  of  .1  bar 
were  held  for  1.5  seconds.  The  ramps  were  continued  until  a  cavita¬ 
tion  event  was  detected.  Twenty  threshold  measurements  were  taken 
at  each  pulse  duration.  PRF.  and  frequency  combination. 


A.  RESULTS  AND  DISCUSSION 


Figures  5.1  and  5.2  are  graphs  of  threshold  pressure  versus  pulse 
duration  taken  at  986  kHz  with  PRFs  of  200  Hz  and  1000  Hz,  respec¬ 
tively.  The  error  bars  in  these  and  subsequent  figures  represent  one 
standard  deviation  of  the  data  above  and  one  standard  deviation  below 
the  mean  threshold  pressure.  It  should  be  noted  that  thresholds  were 
undetectable  at  the  2  (isec  pulse  width  for  the  200  Hz  PRF.  most  likely 
due  to  lack  of  scattered  signal  from  a  two-cycle  burst  at  such  a  low 
PRF.  At  1000  Hz  PRF,  the  thresholds  are  essentially  constant  down  to 
the  10  psec  pulse  duration.  At  5  and  2  psec  pulse  durations,  the 
thresholds  show  a  marked  increase.  The  data  in  Figures  5. 1  and  5.2 
show  good  agreement  with  Reference  6. 
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Figure  5. 1 

Threshold  Pressure  vs.  Pulse  Duration  986  kHs,  200  Hr  PRP 
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Figure  5.2 

Threshold  Pressure  vs.  Pulse  Duration  086  kHz,  1000  Hz  PRF 


Figures  5.3  and  5.4  are  graphs  of  threshold  pressure  versus  pulse 
duration  at  2.17  MHz  for  PRFs  of  200  and  1000  Hz.  respectively. 
There  is  an  essentially  constant  threshold  at  200  Hz  PRF  for  pulse 
durations  of  15  nsec  and  larger,  while  shorter  pulse  durations  indicate 
a  marked  increase  in  threshold  pressure.  In  contrast,  however,  at 
1000  Hz  PRF.  the  thresholds  are  essentially  constant  for  all  pulse 
durations  down  to  2  nsec,  with  a  marked  increase  evident  at  a  pulse 
duration  of  1.5  nsec. 

Figure  5.5  is  a  comparison  of  the  thresholds  obtained  for  the  two 
source  frequencies  at  200  Hz  PRF;  Figure  5.6  is  the  same  at  1000  Hz 
PRF.  These  results  support  theoretical  predictions  by  Flynn  IRef.  41 
and  Apfel  (Ref.  31  that  the  threshold  is  independent  of  frequency  for 
sufficiently  small  cavitation  nuclei.  For  parameters  used  in  this 
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Figure  5.3 

Threshold  Pressure  vs.  Pulse  Durmtlon  2.17  Bias.  200  Hz  PRF 
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Figure  5.4 

Threshold  Pressure  vs.  Pulse  Durmtlon  2.17  Bills,  1000  Hz  PRF 
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Pulse  Duration  (microseconds) 


Figure  5.6 

Threshold  Pressure  vs.  Pulse  Durmtion 
986  kHz  end  9.17  MHz  st  1000  Hz  PRF 


experiment,  threshold  pressure  values  are  essentially  equal  to  Blake 
pressure  thresholds,  which  have  no  frequency  dependence  (see 
equation  2,11).  This  conclusion  is  valud  due  to  the  close  relationship 
between  transient  and  Blake  thresholds  for  the  frequencies  under 
consideration  (see  Figure  2.3). 

Figure  5.7  is  a  histogram  of  pressure  thresholds  (rounded  to  the 
nearest  bar)  measured  for  each  individual  transient  cavitation  event. 
Again  from  equation  (2.11)  above,  there  is  a  direct  relationship 
between  bubble  size  and  Blake  threshold  pressure.  Figures  5.8 
through  5.11  are  histograms  of  threshold  measurements  for  individual 
frequency/ PRF  combinations.  It  should  be  noted  that  the  graphs  may 
be  slightly  skewed  due  to  the  starting  pressure  amplitudes  of  most 
ramps.  Since  our  experience  indicated  that  expected  threshold 
pressures  were  in  the  6-9  bar  range,  source  pressure  ramps  were 
started  in  the  3-4  bar  range.  An  occasional  low  threshold  event  due  to 
the  presence  of  an  unusually  large  nucleus  may  have  been  missed. 
Keeping  this  in  mind,  it  appears  that  all  histograms  except  Figure  5,8 
exhibit  a  roughly  Gaussian  distribution  with  mean  values  in  the  7-8  bar 
range.  This  corresponds  to  bubble  radii  in  the  range  of  0.075-0.09 
|im. 

As  stated  above,  when  larger  pulse  width  measurements  are 
compared  in  Figures  5.5  and  5.6,  the  thresholds  are  independent  of 
pulse  duration.  As  pulse  durations  decrease,  there  is  an  increase  in 
threshold.  This  observation  would  indicate  that  the  larger  pulse 
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Figure  5.7 


Histogram  of  All  Threshold  Pressure  Measurements 
(rounded  to  nearest  bar) 
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Figure  5.8 


Histogram  of  All  Threshold  Pressure  Measurements 
088  kHz.  200  Hz  PRF 
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Figure  5,9 

Histogram  of  All  Threshold  Pressure  Measurements 
986  kHz*  1000  Hz  PRF 
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Figure  5.10 

Histogram  of  All  Threshold  Pressure  Measurements 
2.17  MHz,  200  Hz  PRF 


durations  allow  time  for  a  growth  mechanism,  such  as  rectified 
diffusion,  to  cause  nuclei  to  evolve  to  a  Blake  or  transient  threshold 
radius.  There  appears  to  be  insufficient  time  for  the  growth 
mechanism  to  take  effect  at  shorter  pulse  durations  It  is  postulated 
that,  during  the  rest  time  between  pulses,  nuclei  dissolve  to  their 
original  equilibrium  sizes. 

Further  examination  of  Figures  5.1  through  5.4  reveals  that  there 
is  a  distinct  relationship  between  the  flat  and  sloped  portions  of  the 
graphs.  Specifically,  Figure  5.1  (low  frequency/ low  PRF)  reveals  a 
continuously  decreasing  threshold  pressure  for  increasing  pulse 
duration.  Although  no  additional  data  was  taken  at  pulse  durations 
greater  than  40  nsec,  there  were  essentially  no  thresholds  recorded  at 
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levels  lower  than  the  threshold  determined  at  that  point,  so  it  is 
anticipated  that  the  graph  flattens  out  beyond  this  point.  An  attempt 
was  made  to  evaluate  this  hypothesis  experimentally,  but  problems 
were  encountered  with  the  986  kHz  source  transducer  that  prevented 
any  further  reliable  data  from  being  obtained  at  that  frequency.  Figure 
5.2  (low  frequency/high  PRF)  shows  relatively  flat  threshold 
measurements  for  pulse  durations  of  10  psec  and  greater,  with 
significantly  higher  thresholds  occurring  at  5  and  2  psec  pulse 
durations.  Figure  5.3  (high  frequency/low  PRF)  shows  flat  threshold 
measurements  for  pulse  durations  of  15-20  psec  and  greater,  and 
increasing  thresholds  for  10,  5.  and  2  psec  pulse  durations.  Finally, 
Figure  5.4  (high  frequency/high  PRF)  shows  essentially  flat  threshold 
measurements  down  to  2  psec  pulse  duration.  This  observation 
prompted  an  additional  set  of  10  threshold  measurements  at  2.17 
MHz.  1000  Hz  PRF,  and  1.5  psec  pulse  duration  (3  cycles  per  pulse)  to 
ascertain  whether  there  was  a  rise  in  threshold  below  2  psec.  The 
mean  value  of  the  ten  readings  was  12.9  bar,  with  standard  deviation  of 
2.90,  which  indicates  a  “knee"  in  this  graph  as  well  in  the  vicinity  of 
the  2  psec  pulse  duration.  There  appears  to  be  a  relationship  between 
frequency,  pulse  duration,  and  PRF  which  allows  a  certain  range  of 
pulses  per  unit  time  to  mark  the  break  point  where  shorter  pulse 
durations  begin  to  affect  cavitation  thresholds.  In  Reference  6  it  was 
postulated  that  threshold  became  independent  of  pulse  duration  when 
pulses  contained  more  than  10  cycles.  Our  measurements  show  the 
“knee"  to  be  located  as  Indicated  in  Table  5. 1: 
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TABLE  5.1 


CYCLES  PER  IflflX  TIMS  FOR  VARIOUS  FREQUENCY* 

PRF.  AND  PULSE  DURATION  COMBINATIONS 


f 

PRF 

Pulse  length 

(“KNEE") 

Cycles/pulse 

Cycles/sec 

(PRF  x 

Cycles/pulse) 

986  kHz 

200  Hz 

40  (isec 

40 

8,000 

986  kHz 

1000  Hz 

10  psec 

10 

10,000 

2.17  MHz 

200  Hz 

15  psec 

33 

6,600 

2.17  MHz 

1000  Hz 

2  psec 

4 

4,000 

Although  the  value  of  PRF  x  (cycles/ pulse)  product  varies  from 
4,000  to  10.000.  this  analysis  indicates  that  the  knee  of  the  graph  is 
predictable  within  a  reasonable  range  of  pulse  widths.  In  addition,  our 
choice  of  pulse  widths  was  not  geared  toward  finding  this  break  point 
of  the  graph.  There  may  be  other  pulse  durations  that  more  accurately 
indicate  the  knee  between  those  examined  which  yield  a  narrower 
range  of  cycles  per  unit  time  values.  To  test  the  hypothesis,  a  series  of 
threshold  measurements  were  taken  at  2.17  MHz  and  500  Hz  PRF.  It 
was  anticipated  that  the  knee  of  the  curve  should  occur  in  the  vicinity 
of  the  10  psec  pulse  duration.  Figure  5.12,  which  illustrates  the 
results  of  these  measurements,  shows  that  the  knee  occurs  as 
predicted.  Although  it  is  not  prudent  to  generalize  when  discussing 
mechanisms  as  complex  as  bubble  dynamics  in  sound  fields,  it  is  felt 
that  this  aspect  of  these  findings  bears  further  investigation. 
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Pulse  Duration  (microseconds) 

Figure  5-12 

Threshold  Pressure  vs.  Pulse  Duration  2.17  MHz,  500  Hz  PRF 
B.  RECOMMENDATIONS 

All  experimental  measurements  were  taken  using  the  passive 
acoustic  scattering  method,  which  proved  to  be  an  excellent  detection 
method  when  pulse  widths  were  large.  The  resulting  signal  scattered 
from  the  bubble  “cloud’*  generated  by  the  cavitation  event  resulted  in  a 
clear,  distinct  indication  of  a  transient  event..  As  pulse  durations  were 
decreased,  however,  the  scattered  signal  was  sometimes  weak  and 
lacking  definition.  It  is  essential  that  sonoluminescence  detection 
and/or  active  acoustic  detection  schemes  be  fully  incorporated  into 
the  procedures  to  minimize  any  systematic  errors  that  may  exist  in 
the  passive  detection  system.  The  optical  detection  system  was 
limited  only  by  unavailability  of  luminol  {a  sonoluminescence 
enhancer)  until  the  very  final  weeks  of  the  project.  By  then,  it  was  too 
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late  to  bring  the  system  to  full  operation.  It  is  anticipated  that  this 
system,  when  fully  operational,  will  allow  the  replacement  of  the 
human  element  of  the  current  detection  scheme  with  the  electronics 
system.  As  gathering  of  cavitation  data  via  manual  observation  of  the 
transient  signals  on  an  oscilloscope  requires  exceptional  vigilance  and 
perseverance  on  the  part  of  the  operators,  fully  automating  the 
detection  mechanism  is  essential  to  the  continued  success  of  this 
apparatus. 

It  is  recommended,  when  the  system  is  fully  operational,  that 
additional  threshold  measurements  be  taken  at  specific  pulse  dura- 
tion/PRF/frequency  combinations  in  order  to  confirm  the  cycles  per 
unit  time  effect  noted  in  the  hypothesis  above.  Extremely  low  or  high 
PRFs  (10  Hz,  8  kHz)  should  indicate  whether  the  hypothesis  has 
merit.  Localization  of  a  specific  number  of  pulses  or  period  of 
insoniftcation  which  marks  the  knee  of  the  threshold  curve  may  lead 
to  a  better  determination  of  the  mechanisms  which  allow  the 
cavitation  events  to  fully  evolve. 


APPENDIX  A 


!0 

20 

30 
40 
50 
60 
70 
80 
90 
100 
1 10 
i  20 
130 
l  40 
;  i  50 
1  SO 
!  70 
180 
1 90 
200 
‘.:i0 
220 
230 
240 
250 
260 
270 
260 
290 
300 
2(0 
220 
1 330 
1240 
;2S0 
;280 
370 
250 
393 

m 

410 

<20 

420 
443 
4S3 
463 
<70 
4  SO 
490 

sao 


APPARATUS  CONTROL  PROGRAM  “CAVDAT- 
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THE  PURPOSE  OF  "CAUDAT"  IS  TG  CREATE  AN  ULTRASONIC  POUND  FIELD 
BY  PROGRAMMING  THE  HP  331  4A  FUNCTION  GENERATOR.  USER-DEFINED 
PULSE  BURSTS  ARE  SENT  TO  A  FOCUSSED  TRANSDUCER  WHICH  IS  THEN 
SWEPT  THROUGH  A  USER-DEFINED  AMPLITUDE  RANGE.  THE  USERS  CAN  THEN 
OeWTE  THEIR  ATTENTION  TO  DETECTION  Gc  CAVITATION  EVENTS  DISPLAYED 
ON  AN  OSCILLOSCOPE.  FURTHER.  WHEN  A  CAVITATION  EVENT  IS  DETECTED, 
AND  THE  COMPUTER  IS  MANUALLY  C'JED,  THE  COMPUTER  WILL  DETERMINE 
TRANSDUCER  SOURCE  LEVEL  USING  CONVERSION  FACTORS  PRE-OETER.UNED 
EXPERIMENTALLY.  THE  CRITICAL  USER  INPUTS  ARE  PULSE  REPETITION 
FREQUENCY  i PRF  >  AND  PULSE  DURATION.  THIS  ALGORITHM  WILL  CALCULATE 
THE  PECUIREO  NUMBER  CF  CYCLES  PER  BURST,  N.  TO  ACHIEVE  THE  DESIRED 
WAVEFORM  SUBJECT  TC  THESE  PARAMETERS, 


. . ••VARIABLE  DEFINITIONS*********** . . 

! 

!  Ampout  *  HP  3314A  OUTPUT  VOLTAGE  AMPLITUDE  IN  MILLIVOLTS. 

1  Aocsiari  *  HP  33UA  STARTING  VOLTAGE  AMPLITUDE  IN  MILLIVOLTS. 

’  Gain  -  PQtf£P  AMPLIFIER  GAIN. 

1  Incr  »  HP  33 * 4A  VOLTAGE  INCREMENT  IN  MILLIVOLTS. 

‘  in ivi  •  SU/Tft  INTERVAL  IN  MILLISECONDS. 

1  N  -  NUMBER  OF  CYCLES  IN  PULSE  BURST i  HP  33 UA  IS  IN  N-CYCLE  MODE. 
‘  K  -  COUNTER  INDEX, 

'  L  «  COUNTER  INDEX. 

i  Nu<**b  steps  -  NUMBER  OF  STEPS  IN  VOLTAGE  AMPLITUDE-  SWEEP. 

1  Oofreq  *  TRANSDUCER  OPERATING  FREQUENCY. 

I  P««HK>  -  PEAK  PRESSURE  AT  CAVITATION  EVENT  IN  BAR. 

'  Prt  •  PULSE  REPETITION  FREQUENCY  IN  HZ. 

'  Pulse  •  PULSE  DURATION  IN  MICROSECONDS. 

!  Re*  *  REMAINDER  OF  THE  INDICATED  DIVISION. 

i 

CLEAR  SCREEN 

clear  ?  i clear  hp-ib  interface 

Gn  KBO  ALL  GOTO  1 *70  'ALLOW  FOR  USER  INTERRUPT  Q>  PROGRAM . 


' ♦•♦••♦•••••VARIABLE  DECLARATIONS************ 

l 

INTEGER  K  X .N.Nuko steps 

REAL  Aftoout  .*««ps tart  .Gein.lncr  .Int  *1  .Oofrec  .Prf  .Pulse  .Ren 
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i  •»*»*»* **••••*•»»«*•« *G  THEN'S  I  ON  LINEAR  nfiR.MYr  *********** ***♦♦**#•*♦« 

5:8  1 

5.30  DIM  :0CO  -  .Peak*'  I3W) 


540  i 

c:t0  i  *«*»»••»*»*«»•»»•#»*••**»*•*■» ADDRESSES** *•♦*♦**♦*♦**«'**»*••♦******♦« 

5£0  ' 

5?3  !  IS  rHc  AQGSESS  OF  Tn£  Hr  33i4A  FUNCTION  GENERATOR. 

580  '  73C  IE  THE  ADDRESS  OF  THE  HP  CCC5A  THINK J£T  PRINTER. 

590  ! 

ECO  < 

£ !0  i . *  IMPORTANT* «•»•»***** **•***»*•♦* . 

ECO  ' 

030  ''  PRIOR 'TO  ANY  TEST  RUN ,  UERIFY  POWER  AMPLIFIER  GAIN.. 

S40  ' 

£50  »* . . .****INP'Jt  PULSE  PARAMETERS****'****************** 

550  ' 


s-a  Clear  screen 

5S0  0 ISP  ’CLEA3E  ENTER  TRANSDUCER  OPERATING  FREQUENCY  IN  KHZ.'i 

590  INPUT  Ocfreq 

700  DI3P  “PLEASE  ENTER  POWER  -MFLIFIEF.  GAIN  COO  OR  LESS).** 

7 i 3  INPUT  Gam  •  OPERATOR  INPUTS  THE  GAIN  FOR  r'JTURt  REFERENCE. 

7C0  IF  GiiiY'OOQ  THEN  GOTO  ?40  v 
730  SOTO  770 

740  CIS*  “THIS  GAIN  MAY  RESULT  IN  TRANSDUCER  DAMAGE . “ 

'53  WAIT  2 

7G0  GOTO  730 

770  QISP  “PLEASE  ENTER  CE3IRED  PULSE  REPETITION  FREQUENCY  -IN 

IM»UT  Pp  f 

*•90  Irv.i-n-Prf  >*1000  -CALCULATE  SW,  TP  INTEB'ML  IN  MILLISECONDS. 
933  EISA  “PLEASE  ENTER  DESIRED  PULSE  DURATION  IN  MICRO SECONDS. *  < 

910  INPUT  ?ul»« 

833  N"<  < Pj; ae«Cp f req  1 030  )  'CALCULATE  NUMBER  OF  CYCLES. 

G30  'ECHO  THE  INPUT  PARAMETERS  AND  AlLOW  USER  *0  EDIT  AS  NECESSARY. 
940  PRINT  “TRANSDUCER  OPERATING  FREQUENCY  IS:' lOpfreqt’KHl.* 
eS8  PRINT  “PULSE  REPETITION  FREQUENCY  IS: ‘ jPrf » “M2. • 

360  PRINT  ‘PULSE  DURATION  IS : * tPulaet “MIC3QSECQN0S; * 
e?0  PRINT  “NUMBER  OF  CYCLES  PER  BURST  I$»\«N 
ESS  PRINT  “POuER  AMPLIFIER  GAIN  ISt'iSam 
S30  PRINT  ‘SU/TR  INTERUAL  IS: * t Intv  1 1 “MILL I SECONDS . * 

300  PRINT  “ARE  THESE  PARAMETERS  CORRECT?  (Y/N)* 

S!3  INPUT  SciH 

930  !c  UPC*<DecU)-“N*  THEN  SOTO  570 
J30  ! 

340  '.**•*•**. **.*..********INPUT  sweep  parameters* . 

358 

sse  clear  screen 

970  oisf  -Please  input  sweep  starting  amplitude  in  millivolts. “* 

930  INPUT  Ano start 

938  CISP  “PLEASE  INPUT  AMPLITUDE  INCREMENT  it*  MILLIVOLTS. 

1030  INPUT  tncr 
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-«  .  •.  jk gK Ul SA  t  AH ^C3 SWlAl  'LAlSkl 


».a «- . ad aa aa as s. a 


0 !  0  01  rP  ,£?L£ASs  IN'°UT  NUKEEF  OF  S"F£  F£R  £W£EC " i 

022  INPUT  Nunostscs 

_'3u  IP  And  start '  30vSi3  THEN  3CTC  !C?0 

040  If  Nufi33teDS*Inrr  3000  THEN  5CTC  1070 

US3  IF  Awpstart  +  I  Nunos -cp3*incr- .*  3200  THEN  GOTO  1270 

050  GOTO  1 ! 3C 

070  01 EP  "ENTRY  OF  THESE  PARAMETERS  MAY  RESULT  IN  TRANSDUCER  DAMAGE." 

0sc  wait  : 

C93  6GTC  350 
120  1 

!)C  i ECHO  THE  SWEEP  Parameters. 

::o  t 

130  PRINT  ’  STARTING  AMPLITUDE  IS:  “ ; restart  t  "MILLIVOLTS. " 

140  CRINT  "AMPLITUDE  INCREMENT  IS:  ‘ 'icrj  “MILLIVOLTS.  * 

!50  PRINT  "NUMBER  OF  STEPS  PER  SLEEP  S : "  iNu.-Cstaps 
ISO  PRINT  “ARE  THESE  PARAMETERS  CORRECT’  <Y/N)."i 
170  INPUT  Das 2* 

ISO  IF  L'PC$‘  Dsc2$  '=  “N '  THEN  SOTO  350 

! 90  h"1  1  initialize  storage  field  counter. 

200  1 

213  . *»•♦»**•♦ . ♦•♦•♦♦AMPLITUDE  SWEEP  I NFQRMAT I ON* •♦•♦♦♦»*•«♦♦♦•*♦**• * 

220  i 

230  CLEAR  SCREEN 

240  OUTPUT  ?! 7i  "FR“  sGpfrea » *K2*  *  !aET  TRANSDUCER  OPERATING  FREQUENCY. 

250  OUTPUT  7 1 7 i “M03"  ‘SET  N- CYCLE  MODE. 

260  OUTPUT  7 1 7t “NH*  jNi "EN  ‘  ‘SET  NUMBER  OF  CYCLES. 

270  OUTPUT  7t?j*T:Mn‘.vi>*MS*  ’SET  5U/TR  INTERVAL. 

:se  WAIT  I  ‘ALLOW  ‘SANiiENTS  TO  SETTLE  OUT. 

290  Amaout«A«pslart 

300  FOP  I-i  TO  Nunoateoa  STEP  |  *  COMMENCE  AMPLITUDE  SWEEP. 

310  Ai*CJns(i  l»A«DOUt 

320  OUTPUT  7l?j-AP'iA«o»«c{i >i*MVO*  -SET  FUNCTION  GENERATOR. 

330  Awpout «A«piftc f  I  )♦  I r>cr  'INCREMENT  COUNTER, 

340  U*17  I. 5  ‘ALLOW  TRANSIENTS  TO  SETTLE  OUT. 

350  NE*T  I  ‘CONTINUE  WITH  SWEEP  UNTIL  COMPLETE. 

360  PRINT  "AMPLITUDE  SWEEP  COMPLETE . * 

370  IF  Opfraq-985  THEN  GOTO  1390 

330  IF  Oofreo-2170  THEN  GOTO  1410 

330  PeaKU‘)*><  .3522*U Anpout-Inc-  ‘*i  .2>>) 

480  GOTO  U20 

413  Pea>  <K  >“t .0403*1 < Anpoui-Incr '•;  .2  > » ) 

!  420  PRINT  "SINCE  NO  EVENT  DETECTED  .  THRESHOLD  IS  GREATER  THAN:  *  tPaaMK  >i  "SAP.  * 
1430  GOTO  IS80  lOUERV  USER  FOR  Fw-PTh£R  INSTRUCTIONS. 

1440  » 

1450  , ‘ **«**♦*♦•**♦*♦♦•♦« ♦♦•♦•CAUI TAT  ION  INFCRHAT ION^*****^*»»***« »***•*♦*•* 

1460  *« 

1470  DIFP  -WAS  A  CAVITATION  EVENT  SETECTEO?*; 

1460  INPUT  OecOl 

1490  IF  UPC*(C-ac3*  »«-N*  Th£F4  GOTO  15SS 
1533  IF  UPC*-’0«c3*5-‘Y*  ThEN  GOTO  !Si3 
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1510 
iSZO 
1530 
i  ~“*0 
I  550 
1560 
1  570 
1530 
1530 
1500 
1510 
!  620 
1 630 
i  54® 
1650 
1 660 
1670 
1630 
1650 

1709 

1710 
;  ?2® 
1730 
l  “40 
1750 
>750 
1770 
1733 
1790 
i  60S 
!3ie 
1 520 
i230 
1340 
1650 
1350 
1370 
1650 
1830 
1900 
>910 
1320 
1330 
1349 
I9S0 
i960 
>979 
i  380 
1399 
2000 


IF  Opfreq-983  ThEN  GOTC  1520 
IF  0prraqcZi7C  THEN  GOTO  ',550 
Fear-  <  K  )»?  .0522  «<  Ar.p  ir.ci  I  )*>  .  Z  ? )  > 

GOTO  1550 

Fear  K  ’  =  .3403* hrc : nci  I  )*•.  .2  ;  >  5 

PRINT  "PEAK  PRESSURE  AT  CAVITATION  EVENT 


NUMBER- jNi ‘IS 


;P«aMK >i“SAh. * 


K=»K  + 1 

CIEF  '00  TuU  >.jI3H  TO  RUN  AGAIN  WITH  THE  SAME  PARAMETERS’  *  t 
INPUT  Dec 4$ 

IF  UFC$<  Dec 4$  >=“  Y  >'  THEN  bOTO  1520 
IF  UFCii  Dec4$  .'c”N“  THEN  GOTO  1740 

OUTPUT  7i 7 1 ‘PR"  ‘RESET  HP  3314A  FUNCTION  GENERATOR. 

WAIT  50  1 ESTABLISH  QUIESCENT  PERIOD  BETUEEN  TRIALS. 

2EEF  .56,1  i CUE  USER  THAT  ANOTHER  TRIAL  IS  "BOUT  TO  START. 
5EEF  1 322 . 08 , . 5  * 

PEEP  ‘546. 21,. 5 
BEEP  >750.35,. ?S 
BEEF  1 332.23 , . 25 
EEEP  1790.36,1 
GOTO  1240 


•  . * . *.,****♦• ,.RJUNT-GUT  OPTION**** . 

i 

CI3F  “ WOULD  VGU  lIkE  a  HARO  cop v  of  PARAMETERS  and  RESULTS’*! 

INPUT  Dec5* 

IF  UPCSi  DacSS  )»*v*'  THEN  GOTO  1789 
IF  UPCJf DacSS  i-*N*  THEN  GOTO  1920 
PRINTER  IS  “02  * IDENTIFY  ADDRESS  OF  TH INKJET  PRINTER. 

PRINT  ‘TRANSDUCER  OPERATING  FREQUENCY  IS:  iOptfrftQt‘Krt:.* 

PRINT  'NUMBER  OF  CYCLES  PEP  BURST  IS:  *  tN 

PRINT  -COWER  AMPLIFIER  "AIN  !S:*tSsin 

“PINT  -SW.-TR  INTERVAL  IS  t  “ » Intvl  t  ‘MILLISECONDS, * 

PRINT  -PULSE  DURATION  IS: * jFulset ‘NICRCSECONQS.* 

PRINT  "PULSE  REPETITION  FREQUENCY  IS: * iPrf  >  “HZ . * 

POR  L-i  TO  a-|i  STEP  I  1  CREATE  DO  LOO®  FOR  PRINTING  OF  RECORDED  E'.SNTS. 
PRINT  -PEAK  PRESSURE  AT  CAVITATION  EVENT  NUMBER*  iU  *i5* *°6d» » L  s  ‘GAR.  * 
NEXT  L  ‘CONTINUE  UNTIL  COMPLETE. 

PRINTER  IS  CRT  I  RETURN  TO  CRT 
1 

•  « . »•***.* *.***CONT TESTING  IF  NECESSARY*. **•*.»•••*«*•••**• 

i  .  . 

OISP  ‘00  YOU  wrSH  TO  ALTER  PARAMETERS  ANC  RUN  AGAIN’*! 

INPUT  Oac3» 

IF  uPCt'Oacb *)•***  THEN  GOTO  I960 
IF  iJPCttOecb*  ‘-*N*  THEN  GOTO  2050 
OU  <  PUT  7l7i  *PR*  IRESET  hP‘33UA  FUNCTION  GENERATOR. 

WAIT  69  'ESTABLISH  QUIESCENT  PERIOD  BETWEEN  TRIALS. 

BEEP  9“S.SG,1  <CUE  USER  THAT  ANOTHER  TRIAL  IS  ABOUT  TO  S'AST. 

BEEP  ! 302.06, .5 
BEEP  1545.22 , .5 
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2310  BEEP  i ”SC . 3S , . ?5 
2C20  BEEP  '302.08  ,.25 
2330  SEEP  ! ?30 . 36  , ! 

2343  c'TC  £?0 

2350  PRINT  "PRCGRHtt  COMPLETE. * 

2050  OUTPUT  7 ! 7 r?R *  'RESET  HP  33 UR  FUNCTION  3ENEEA7CR. 
2373'  £N0 
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APPENDIX  B 


PVDF  HYDROPHONE  CONSTRUCTION  AND 

fiEBBRA L  INFORMATION 

A.  CONSTRUCTION  DETAILS 

Our  desire  for  small  size  led  us  to  choose  a  large  plastic  syringe 
and  a  large-diameter  hypodermic  needle  as  the  hydrophone  housing. 
This  gave  us  a  convenient  method  for  insulating  the  “hot"  lead  from 
the  grounded  elements,  as  well  as  providing  ample  storage  area  close 
to  the  active  elements  for  the  preamplifier  circuitry.  Our  primary  dif¬ 
ficulty  in  construction  was  working  with  small  sizes  of  PVDF  film. 
After  numerous  unsuccessful  attempts  to  cut  the  PVDF  to  a  shape 
which  allowed  for  good  electrical  connection  while  minimizing  active 
element  size  (1  mm2),  we  moved  to  a  PVDF  strip  approximately  2  mm 
wide  and  about  4  cm  In  length.  We  reduced  the  bottom  (positive)  side 
of  the  film  by  etching  away  the  PVDF  using  a  Q-tip  soaked  in  ferric 
chloride,  which  gave  us  a  marginally  accurate  method  of  reducing 
element  size  while  keeping  a  capability  to  have  a  reliable  electrical 
connection.  The  top  (ground)  was  not  etched.  The  PVDF  strip  was 
then  centered  on  a  3.6  mm  x  6.4  mm  brass  plate  attached  to  the  "hot" 
lead.  The  barrel  of  the  needle  was  insulated  with  several  layers  of  heat 
shrink  tubing  to  prevent  any  of  the  “hot"  side  from  short-circuiting  to 
ground.  The  remaining  lengths  of  the  PVDF  strip  were  attached  to  the 
insulation  along  the  needle  barrel  and  connected  to  the  surface  of  the 
needle  with  nickel  paint.  Gaps  between  the  PVDF  strip  and  needle 


were  filled  with  epoxy  adhesive  to  add  stability.  Resistance  measure¬ 
ments  were  taken  to  ensure  continuity  of  connections  and  that  no 
short  circuits  existed.  q'he  unity  gain  preamplifier  (voltage  follower) 
was  constructed  and  installed  in  the  syringe  and  electrical  connec¬ 
tions  made  to  the  hot  lead  and  external  ground  connections.  The 
entire  apparatus  was  insulated  by  being  coated  with  “Plasti-Dip."  See 
Figures  B.i  and  B.2. 

B.  SENSITIVITY  CALCULATION  [Ref.  28] 

E  =  Vo/t  where  E  =  Electric  Field  developed 

Vo  =  Voltage  Induced 
t  =  Thickness  of  Piezoelectric  film 

Also,  E  a  g3aX3  where  g33  is  a  constant  indicating  field/applied  stress 

X3  =  Applied  mechanical  stress 

Mo  =  Pressure  Sensitivity  =  V/P  ^3^3t  =  g33t 

Numerically, 

M0  *  g33t  =  339  x  10-3  x  (9  x  10-«  m) 

=  3.U51  x  IQ"6  =  3.051  x 

M0  =  20  log  10  (3.051  x  10-6)  dB  re  lV/Pa  =  -110.31  uB  re  lV/Pa 
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thrift 
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Exploded  Drawing  of  FVDF  Hydrophone 


This  calculation  indicates  that  transducer  pressure  sensitivity  is 
frequency  independent.  It  should  be  mentioned  that  the  calculated 
sensitivity  is  contingent  upon  the  assumption  that  the  nickel  paint 
used  to  mount  the  PVDF  to  the  brass  mounting  block  effectively 
clamps  the  active  face  in  the  (1)  and  (2)  directions,  thus  confining  the 
face  to  operation  solely  in  the  thickness-expander  mode.  See  Figure 
B.3. 

-eo 
15 

a. 

>  -90 

ci 


3  “100 


c 

<1> 


-120 

0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0 

Frequency  (MHz) 

Figure  B.3 

PVDF  Hydrophone  Sensitivity  Calibration  Curves 

C.  DESCRIPTION  OF  TEST  AND  CALIBRATION  TECHNIQUE 

We  chose  the  reciprocity  technique  for  calibration  of  our 
hydrophone  sensitivity  because  calibrated  hydrophones  in  the 
megahertz  range  were  not  available  for  comparison.  Our  frequency 
band  of  interest  was  0.5  MHz  to  4.0  MHz.  For  the  reciprocity 
calibration,  we  used  two  Panametric  V323  transducers  as  the  trans- 
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mitting  and  reversible  elements  of  the  experimental  setup.  These 
transducers  are  designed  for  operation  at  2.25  MHz  and  provided 
smooth  response  through  the  frequency  range.  The  calibration  set  up 
was  as  noted  in  the  equipment  diagram  included  as  Figure  4.1.  The 
distance  between  the  transducers  was  30  cm,  which  is  well  into  the 
far  field  for  all  frequencies  of  interest.  Far-fleld  calculation  for  the 
maximum  frequency  is  as  follows: 

r  =  [Ref.  25] 


where  L  =  5  x  10-3,  m  =  largest  extent  of  source  active  element,  X  = 
wavelength  at  4.0  MHz  in  fresh  water,  r  =  minimum  extent  of  far  field, 
and 

(5  x  10-3m)2  A  . flon  ,.on 
r  =  4(3.70  x  lO^m)  =  m  =  16-89  cm 


We  took  Voltage  Level  (in  dBV)  and  input  current  (in  milli- 
amperes)  readings  from  the  HP  3585A  Spectrum  Analyzer  and  HP 
3478A  Digital  Multimeter  in  accordance  with  standard  reciprocity 
techniques. 

Open-circuit  voltage  sensitivity  of  our  PVDF  hydrophone  was 
determined  using  the  following  equation: 


M0x* 


(2Xt\ 

tpocj 


VxVx 


Y\l/2 


vr!r  ) 


[Ref.  251 


When  converted  to  sensitivity  level,  this  calculation  became: 
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D.  DIRECTIVITY  EFFECTS/EXCITATION  OF  FIRST  THICKNESS 

MODE 

1.  Analysts  of  Directivity  Effects 

One  possible  cause  of  the  difference  between  actual  and  pre¬ 
dicted  sensitivities  may  be  the  change  in  the  directivity  of  the 
hydrophone  as  frequency  increases.  As  the  transmitting  transducer 
frequency  is  increased,  both  the  hydrophone  and  transducer  become 
more  directional,  implying  more  efficient  transmission  by  the  trans¬ 
ducer  and  more  efficient  reception  by  the  hydrophone  of  the  acoustic 
power  in  a  given  direction.  The  trade-off  is  that  the  main  lobes  of 
both  the  transducer  and  the  hydrophone  become  narrower.  Using  the 
approximation  for  a  piston-like  source  {Ref.  251, 

D  »  1/4  k2LiL2  where  Li  s  2  mm  and  L2  s  5  mm 
k  =  (o/c  =  2rcf/c 

At  500  kHz.  DI  =  10  log  D  -  10.5  dB 

These  “corrected"  calculated  values  are  plotted  as  such  in  Figure  B.3. 

2.  Calculation  of  First  Thickness  Mode  of  Brass  Mounting 

Elitfi 

It  was  originally  believed  that  the  brass  mounting  plate  would 
behave  as  if  it  were  perfectly  rigid  in  the  frequency  range  of  interest. 
Measurements  led  us  to  believe  that  this  may  not  be  the  case.  For  this 


reason,  it  was  decided  to  determine  if  the  first  thickness  mode  of  the 
brass  mounting  plate  was  being  excited  somewhere  in  the  range 
between  0.5  MHz  and  4.0  MHz.  The  first  thickness  mode  will  be 
excited  with  t  =  X/2  or  when  X  =  2t.  For  a  thickness  of  1.3  mm,  X  = 
2.6  mm.  Since  c  =  fX,  f  =  c/X: 

f  =  c/X  =  (4700  m/s)/2.6  x.10-3  m  s  1.81  MHz 

It  becomes  evident  that  it  is  likely  that  the  first  thickness 
mode  is  being  excited  to  some  extent. 

E.  DISCUSSION 

As  noted  above,  theoretical  sensitivity  was  calculated  assuming 
that  the  longitudinal  and  transverse  modes  were  clamped  and  that 
only  the  thickness  mode  was  excited.  Data  showed  little  agreement 
with  these  original  theoretical  calculations,  however.  The  PVDF 
hydrophone's  active  face  was  considerably  larger  than  the  wavelengths 
of  the  frequencies  of  interest,  which  makes  the  receiving  beam  pat- 
tern  i.^ore  narrow  as  the  frequency  increases.  This  directivity  also 
adds  considerably  to  the  sensitivity  when  sound  is  received  “on  axis." 
Directivity  and  Directivity  Index  as  calculated  previously  in  this  paper 
contributed  greatly  to  the  sensitivity  of  the  hydrophone.  As  can  be 
seen  from  the  graph  of  Pressure  Sensitivity  vs.  Frequency  (See  Figure 
B.3),  there  is  good  agreement  between  theoretical  and  experimental 
values  at  the  lower  end  of  the  frequency  range  (up  to  1.5  MHz).  At 
higher  frequencies,  additional  variables  are  introduced  which  limit 
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agreement.  A  likely  contributor  to  the  higher  frequency  disparity 
might  be  slight  misalignment  of  the  hydrophone.  As  frequency  goes 
up.  the  beam  pattern  becomes  increasingly  narrow.  Slight  misalign¬ 
ment  could  result  in  poorer  reception  of  the  transmitted  acoustic 
energy.  Lastly,  as  alluded  to  previously,  it  is  possible  that  the  first 
thickness  mode  of  the  brass  mounting  plate  is  being  excited  in  the 
frequency  range  of  interest. 


89 


LIST  OF  REFERENCES 


1.  Ross.  D..  Mechanics  of  Underwater  Noise,  p.  202.  Pergamon 
Press.  Inc.,  1976. 

2.  Strasberg,  M.,  “Onset  of  Ultrasonic  Cavitation  Tap  Water,"  Journal 
Of.thfijteflUgtteal  SflCtetE...QLAmsrtea.  V.  31,  pp.  1 63-176,  February 
1959. 

3.  Apfel,  R.  E.,  “Possibilitv  of  Microcavitation  from  Diagnostic  Ultra 
sound,"  IEEE  Transactions  on  Ultrasonics.  Ferroelectrics,  and 
Frequency  Control,  v.  UFFC-33,  pp.  140-146,  March  1986. 

4.  Flynn,  H.  G„  “Generation  of  Transient  Cavities  in  Liquids  by 
Microsecond  Pulses  of  Ultrasound."  Journal  of  the  Acoustical 
Society  of  America,  v.  72,  pp.  1926-1932,  December  1982. 

5.  Roy,  R.  A..  Atchley,  A.  A..  Crum,  L.  A.,  Fowlkes,  J.  B.,  and  Reidy. 
J.  J.,  “A  Precise  Technique  for  the  Measurement  of  Acoustic  Cavi¬ 
tation  Thresholds  and  Some  Preliminary  Results."  Journal  of  the 
ACOUStlC.aL.Sacle.ty-Ql,,,^SIl£a.  V.  78.  pp.  1 799-1 805.  November 
1985. 

6.  Atchley.  A.  A.,  Frizzell.  L.  A.,  and  Apfel.  R.  E..  “Thresholds  for 
Cavitation  Produced  in  Water  by  Pulsed  Ultrasound."  submitted  for 
publication  in  Ultrasonics. 

7.  Walton.  A.  J..  and  Reynolds.  G.  T.,  “Sonoluminescence,"  Advances 
in  Physics,  v.  33.  pp.  595-660.  1984. 

8.  Neppiras.  E.  A.,  and  Noltingk,  B,  E.,  "The  Onset  of  Cavitation  Pro¬ 
duced  by  Ultrasonics,"  Proceedings  of  the  Physical  Society  of 
London  Section  B.  v,  64.  pp.  1032-1038,  1951. 

9.  Harvard  University  Acoustics  Research  Laboratory  Technical 
Memorandum  No.  12.  The  Onset  ofCavltatlQnJn.  Liquids;!,  by  F.  G. 
Blake.  Jr.,  p.  28,  2  September  1949. 

10.  Rayleigh,  Lord.  “On  the  Pressure  Developed  in  a  Liquid  During  the 
Collapse  of  a  Spherical  Cavity."  Philosophy  Magazine,  v.  34,  pp. 
94-98,  August  1917. 

11.  Plesset,  M.  S..  The  Dynamics  of  Cavitation  Bubbles,"  Journal  of 
Applied  Mechanics,  v.  16,  pp.  277-282,  September  1949. 


90 


12.  Poritsky,  H.,  “The  Collapse  or  Growth  of  a  Spherical  Bubble  or 
Cavity  In  a  Viscous  Fluid,"  Proceedings  of  the  First  U.S.  National 
Congress  on  Applied  Mechanics,  ed.  E.  Sternberg  (New  York),  pp. 
813-821,  1952. 

13.  Noltingk,  B.  E..  and  Neppiras,  E.  A.,  “Cavitation  Produced  by 
Ultrasonics,”  Proceedings  of  the  Physical  Society  of  London  Sec¬ 
tion  B.  v.  63.  pp.  674-685,  1950. 

14.  Lauterborn,  W.,  “Numerical  Investigation  of  Nonlinear  Oscillations 
of  Gas  Bubbles  in  Liquids,"  Journal  of  the  Acoustical  Society  of 
America,  v.  59,  pp.  283-293.  February  1976. 

15.  Crum,  L.  A..  “Rectified  Diffusion,"  Ultrasonics,  v.  22.  pp.  215-223, 
September  1984. 

16.  Harvey,  E.  N.,  Barnes,  D.  K.,  McElroy,  W.  D.,  Whitely,  A.  H.,  Pease. 
D.  C..  and  Cooper.  K.  W.,  “Bubble  Formation  in  Animals."  Journal 
of  Cellular  and  .Conmaratlve.PhyslQlQgy,  v.  24.  pp.  1-40,  1944. 

17.  Grifflng,  V..  “Theoretical  Explanation  of  the  Chemical  Effects  of 
ultrasonics."  JamnaLaf.i^ismlcal.,  Physics,  v.  18.  pp.  997-998. 
1950. 

18.  Sehgal,  C.,  Steer.  R.  P.,  Sutherland,  R.  G..  and  Verrall,  R.  E.. 
“Sonoluminescence  of  Aqueous  Solutions,"  Journal,  of  Physical 
Chemistry,  v.  81.  pp.  2618-2620.  December  1977. 

19.  Sehgal.  C„  Sutherland.  R.  G.,  and  Verrall.  R.  E..  “Optical  Spectra 
of  Sonoluminescence  from  Transient  and  Stable  Cavitation  in 
Water  Saturated  with  Various  Gases."  Journal  of  Physical 
Chemistry,  v.  84.  pp.  388-395,  February  1980. 

20.  Sehgal,  C„  and  Verrall,  R.  E„  "A  Review  of  the  Electrical  Hypo¬ 
thesis  of  Sonoluminescence,"  Ultrasonics,  v.  20.  pp.  37-39. 
January  1982. 

21.  Vaughan,  P„  and  Leeman,  S.,  "Some  Comments  on  Mechanisms  of 
Sonoluminescence."  Acustlca.  v.  59,  pp.  279-281.  1986. 

22.  Crum,  L.  A.,  and  Reynolds.  G.  T.,  "Sonoluminescence  Produced  by 
Stable  Cavitation."  Journal  of  the  Acoustical  Society  of  America,  v. 
78.  pp.  137-139,  July  1985. 

23.  Daniels.  F.,  and  Alberty.  R.  A.,  Physical  Chemistry.  4th  ed..  p.  328. 
John  Wiley  &  Sons.  Inc.,  1975. 


91 


24.  Clay,  C.  S.,  and  Medwin,  H.,  Acoustical  Oceanography:  Principles 
and  Applications,  pp.  184-204,  John  Wiley  &  Sons,  Inc.,  1977. 

25.  Kinsler,  L.  E„  and  others.  Fundamentals  of  Acoustics.  3rd  ed..  pp. 
389-391,  John  Wiley  &  Sons,  Inc..  1982. 

26.  Dunn,  F.,  Averbuch,  A.  J.,  and  O’Brien,  W.  D.,  “A  Primary  Method 
for  the  Determination  of  Ultrasonic  Intensity  With  the  Elastic 
Sphere  Radiometer,"  Acustlca.  v.  38.  pp.  58-61,  1977. 

27.  Lindsay.  R.  B..  MeshanteaLRadlatton.  pp.  217,  292-293.  McGraw- 
Hill  Book  Co..  1960. 

28.  Pennwalt  Corporation.  KYNAR™  Piezo  Film  Technical  Manual, 
1983. 


92 


INITIAL  DISTRIBUTION  LIST 


Nfti..  Copies 

1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library  Code  0142  2 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

3.  Dr.  L.  E.  Hargrove  1 

Office  of  Naval  Research 

Physics  Division  -  Code  1112 
80u  N.  Quincy  Street 
Arlington,  VA  22217 

4.  Commander,  Naval  Coastal  Systems  Center  1 

Dr.  Elan  Moritz  -  Code  4120 

Panama  City,  FL  32407-5000 

5.  Professor  A.  A.  Atchley,  Code  61  Ay  4 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

6.  Professor  S.  R,  Baker,  Code  61Ba  1 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

7.  Professor  S,  L,  Garrett,  Code  61Gx  1 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

8.  Chairman,  Department  of  Physics,  Code  61  1 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

9.  LCDR  R.  L.  Bruce  3 

6002  Crocus  Court 

Alexandria,  VA  22310 

10.  LT  R.  D.  Middleton.  Jr.  5 

141  John  Street 

Hauppauge,  NY  11787 


93 


1 1 .  Amy  and  Aaron  Bruce 
213  Brookdale  Drive 
League  City,  TX  77573 

12.  CAPT  Robert  Perron 
Department  of  Physics 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

13.  Mr.  and  Mrs.  R.  D.  Middleton,  Sr. 
141  John  Street 

Hauppauge,  NY  11787 

14.  Maris  R,  Zanger 

300  Glenwooa  Circle  #164 
Monterey,  CA  93940 


94 


